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PREFACE 

This book is a republication of a paper in three parts, of which the 
first part was written by me and the second and third parts by me 
in conjunction with Miss Hilda P. Hudson. They were published 
by the Royal Society in its ‘ Proceedings/ A, 1916 and 1917. 
Since that time, I have not been able to continue these studies. 

The whole of the subject was also dealt with by the Finite 
Calculus at the end of my book “Prevention of Malaria/' 2nd 

, ' «^ ) but this Addendum was written 

in such a hurry and under the press of so much other work, that 
it is not very clear. 

The earlier chapters of the same book, especially Sections 26, 

2/ and 28, also contained an elementary kind of pathometrv for 

non-mathematical readers, which does not pretend to be exact 

ut of which copies can be obtained from my Personal Secretary 
at this Institute. J 

I should add here that the studies contained in our original 
paper published by the Royal Society have been followed up by 
Prof. YY. 0. McCormack and Colonel A. J. McKendrick in another 

mTiSt Th^ 006 !? 1 ? 1 ? W R ° yaI S ° ciety ’’ A ’ Vol ‘ 115 > 

p“ y Yor5£, b / an ^, S* P “ 4 Science 

Meuxan° u s „ r diheteroxenic parasitism is not dealt with in this 
book, but was briefly considered in my “ Prevention of Malaria ,” 2nd 

A VT’tt e6 °. m° and U) The anal yt>cal problem was solved bv 

The s,/ “t NatUI6 ’ • 8 “* Feb - 1912 ’ P- 497 > “d ^where. ' 

L W. Hack// • C0 “P lcat \ d b 5' the fact recently discovered by 

the carriers. be ^ Ween — of malaria and the n“ 
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Section [ 

Introduction.—It is somewhat surprising that so little mathematical 
work should have been, done on the subject of epidemics, and. indeed, 
on the distribution of diseases in general. Not onlv is the theme 
of immediate importance to humanity, but it is one which is funda¬ 
mentally connected with numbers, while vast masses of statistics 
ha\e long been awaiting proper examination. But, more than this, 
many and indeed the principal problems of epidemiology on which 
preventive measures largely depend, such as the rate of infection, 
the frequency of outbreaks, and the loss of immunitv. can scarcely 
ever be resolved by any other methods than those of analysis. For 
example, infectious diseases may perhaps be classified in three groups: 
(1) diseases such as leprosy, tuberculosis, and (1) cancer, which 
fluctuate comparatively little from month to month, though they 
may slowly increase or decrease in the course of years : (2) diseases 
such as measles, scarlatina, malaria, and dysentery, which, though 
constantly present m many countries, flare up in epidemics at frequent 
intervals ; and (3) diseases such as plague or cholera, which disappear 
entirely after periods of acute epidemicity. 

To what are these differences due '? Why, indeed, should epi¬ 
demics occur at all, and why should not all infectious diseases belong 
to the first group and always remain at an almost flat rate '? 



these phenomena there must be causes which are of profound 


importance to mankind and which probably can be ascertained only 
by those principles of careful computation which have yielded such 
brilliant results in astronomy, physics, and mechanics. Are the 
epic enncs in the second class of diseases due (1) to a sudden and 
simultaneous increase of infectivity in the causative agents living 
in a ected persons , or (2) to changes of environment which favour 
t eir dissemination from person to person ; or (3) merely to the 
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increase of susceptible material in a locality due to the maduaI loss 
ot acquired immunity in the population t lie re ; or to similar or 
other causes { And why should diseases of the third class disappear, 
as they undoubtedly do. and diseases of the first class remain so 

♦ v 

persistently t —all questions which immediately and obvious]v 
present themselves for examination. 

The whole subject is capable of study by two distinct met hoc Is 
which an* used in other branches of science, which arc complementary 

of each other, and which should converge towards the same results 

* 

— the a posteriori and the a priori methods. In the former we 
commence with observed statistics, endeavour to fit analytical laws 
to them, and so work backwards to the underlying cause (as done 

t i ' 

in much statistical work of the day) : and in the latter we assume a 
knowledge of the causes, construct our differential equations on 
that supposition, follow up the logical consequences, and finally 
test the calculated results by comparing them with tin 1 observed 
statistics. 

Apparently the first (( posteriori work of importance on epidemics 
was a communication by T)r. AVilliam Farr in I SCO (1). in v 
maintained, in connection with an epidemic of cattle plague, that 
the course of the epidemic would follow a curve ot which the third 
difference of the logarithm was a negative constant. In ]SD>--74. 
Dr. (J. H. Evans (2) endeavoured to extend this law to other epidemic- 
but ended hv expressing disappointment with his results. Abac 
recently, however. Dr. J. Brownlee has continued the work in a sern*> 

X * 

of papers (3-8). In the first of those ho said that Dr. Farrs 
curve is the normal curve of probability, but. alter studying a 
number of epidemics, found that Pearson s Type IV fitted better, 
and he concluded (p. 507 ) that " an epidemic is an organic phenome¬ 
non, the course of which seems to depend on the acquisition by an 
organism of a high grade of nifectivity at the point where the 
epidemic starts, this infec-tivity being lost from that period till the 
end of the epidemic at a rate approaching to the terms of a geo- 

■ogression." He admitted, however (p. 317). that " other 
factors, which are not clear, seem to come into play : but added 
“ that the epidemic ends because of the lack of susceptible person- 
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Firs no evidence in its favour, either from the form of the curve or 
from the facts " (see also pp. 500 and .501). It is obvious 
from the mere examination of many curves of epidemics that they 
are often remarkably symmetrical bell-shaped curves, which, 
however, frequently tend to fall somewhat more slowly than tliev 
rose : and in a later paper (4) Brownlee emphasises this feature of 
symmetry, and says (p. 2) that " the deduction from this phenome¬ 
non is direct and complete, namely, that the want of persons liable 
to infection is not the cause of the decay of the epidemic. On no 
law of infection which I have been able to devise would such a cause 
permit of epidemic symmetry. The fall must in all such cases lie 
much more rapid than the rise, though, on the contrary, when asym¬ 
metry is markedly present the opposite holds." In still later 

papers (6. 7) he gives much evidence to show that measles and 
smallpox have distinct periods of recrudescence. 

So far as I can ascertain, a priori researches on epidemiology were 
first commenced by myself since 1899 in connection with malaria. 
In 1904 I read a paper (9) concerned with the random migration of 
mosquitoes, a subject of vital importance in flu* tlieorv of malaria 
which was subsequently dealt with, at my suggestion, by Prof. K 
Pearson, whose researches were then employed by Brownlee (•>). 


Subsequently I endeavoured to find the malaria equations by a 
prion reasoning (10) ; and in the second edition of mv book (11). 
extended this method to a preliminary general ” Theory of Happen¬ 
ings, employing chiefly the Finite Calculus (which is useful for 
malaria). Quite recently I published a very brief note (13), recordim- 

further advances ; and Brownlee added another (7), showino- how 
he arrived at some of his results. * 

The present paper is a much more advanced and general develon- 
ment of the work commenced in Section 66 (14)° of the second 
edition of my book (11, p. 676). It deals merely with the curves 
w nc ai e theoretically obtained when we suppose that the infecti vitv- 
ratio remains constant or proportional to the number of individuals 
a ea y a ectec, while, simultaneously, some of these are constantly 
re\er ing, or losing immunity, and while both the affected and 
the unaffected groups are subject to special rates of birth, death, 
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immigration and emigration—as defined in the following section. 
This thesis does not consider progressive decline of the infectivity ; 
yet we shall see in Section YII that, contrary to jwhat Dr. Brownlee 
supposed, it does often yield curves which are bell-shaped, nearly 

symmetrical, or. when otherwise, decline more slowly than they rose. 

* ' «/ ^ 

that is. yields curves of the type frequently found in epidemics and 
does so without demanding any other cause for the decline than the 

C v 

exhaustion of susceptible material. 

The thesis is stated in terms which are certainly academical, but 

into which it has been resolved only after very careful thought. It 

has not only been dealt with in considerable detail, but some of the 

% 

more immediate consequences have been followed out—because I 
think that it must first be examined and used as a standard before we 
can proceed to discuss variations due to changes in the infectivity. 
As will Ire seen, the mathematics presents few difficulties ; but. 
owing to the nature of the subject, I have thought it wise to labour 
the proofs somewhat more than some readers will require. 

Beyond this point I have not endeavoured to go, because the appli¬ 
cation of the equations to numbers of known epidemics (which will 
ultimately Ire required) can be made only at great length. The 
paper is therefore pur ely theoretical; but it is one which I think 

will Ire needed for future studies. 

Dr. Farr’s communication (1) is not republished in his collected 
works (‘ Vital Statistics,’ edited by N. A. Humphreys, London, 
1885), but has recently been recovered with some difficulty by Dr. 
Brownlee. It is interesting to note that he seems to have attributed 
the decline of epidemics partly to the attenuation of infectivity, 
but also partly to the fact that the individuals left are less suscep¬ 
tible of attack, either by the constitution or hygienic conditions, 
than those destroyed.” The descriptions of liis methods furnished 
by himself and by Evans are, however, almost unintelligible ; but 
it is evident from the calculated series by law which he gives, 
that the third difference (not the second, as Brownlee said) of the 


logarithms of the series is a negative constant. 

Some of Brownlee’s equations are also difficult to interpret in 
terms of the numbers of individuals affected. I have therefore 
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0 


thought it best to proceed in this paper entirelv on the basis of mv 

■I 

own previous work. 


Section II 




mu before us is as f 



f>ws. Suppose tliat we have a 
population of living tilings numbering V individuals, of whom a 

number Z are affected by someth inu (such as a disease), and the re¬ 
mainder A are not so affected: suppose that a proportion h . dt 
of the noil-affected become affected in every element of time dt. 


and that, comeisclv. a proportion r . dt of the affected become 
unaffected, that is, ie\ert m evorv (dement of time to the non- 


affected group : and lastly, suppose that both groups, the affected 
and the non-affected. are subject also to possibly different birth-rates 
death-rates, and immigration and emigration rates in an element of 


time ; then what will be the number of affected individuals, of new 


cases. and of the total population living at anv time / { 

For the solution of this and the subsidiary problems T have ventured 
to suggest the name ** Theory of Happenings.** ft covers main 


cases which occur not only in pathometrv but in the analysis of* 

_ j • . , * » 


11 


questions connected with statistics, demography, public health, 
the theory of evolution, and even conuneree. politics, and statesman¬ 
ship. The name pathometry (pathos, a happening) was previously 
suggested by myself in antithesis to nosometry (tiosos. a disease) 
for the quantitative study of parasitic invasions in the individual. 




(i) Let ndt, nidi, idt, edt denote respectively the nativitv. mortalitv. 

_ 1 * . * 1 i * 


immigration, and emigration rates of the non-affected part of the 
population in the element of time dt : and N dt, M dt, Idt. Edt denote 
the similar rates among the affected part. Then, as argued in mv 


previous writings and as will be easily seen, the problem before us 
may )e put in the form of the following system of differential 

equations :— 

p = vft -t v 2. 


d P = ( n —m -j- i 



; • A + (N — M 4- I 


<? 





(n 



/ 

* I j , 


m 


hit . A h ( 


e h)dt • A (A -j- r)r// . Z 


( 1 ) 


At 




kA 




r 


I - E - 


kA 


-- r 





(3) 


*^A 


t/V 
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Here <1V consists only of the, variation-elements n, m, i. e. N. M. I. E. 
with their proper signs, while <IA and dZ contain also the happetuntj- 
vlnnent h and the re version-element r. Obviously dV — dA d T Z. as 
it should do. the elements h and r disappearing in the summation 
because they denote only change of condition and not of beimr. 
I he number of the non-affected who become affected in the element 
ot time is hiU . A, and. of course 4 , this number pass over from group 
A to mo up Z—as shown in equations 2 and 3. Similarly rdt . Z 
is the number of the affected who revert in the element of time and 
therefore pass over from group Z to group A. The equations are 
ml quite s ymmetrica l, since N. which should appear in equation 3. 
appears in equation 2—because the progeny of the affected group, 
namely, XVfr. Z. Avill generally be born not affected—an important 
fact which modifies almost all the results. 

(ii) The variation-elements h. m, i, e, and N, 51. I, E, may some¬ 
times be functions of time, especially if the considered events extend 
over long periods : but it will quite suffice at present to take them 
as being constants. If v r e have sufficient data regarding them* we 
can generally calculate them bv the methods of Section VIII for 
whatever small unit of time we adopt : otherwise they must be 
conjectured or assumed to begin with, and be then ascertained by a 
comparison of the integrated equations with known facts. It is 
convenient to write 

r EE it — m + i — e, 

V = N — + I - E. (4) 

Tn some kinds of happening which have no marked effect on the 
birth-rates, death-rates, and immigrations such for instance as 
mild maladies, enlistments in civil professions, conversions to religious 
or political parties, entries into unlimited societies or trades-unions, 
etc.—-wo mav have n — N. m = M. i = I, e = E. and consequently. 


H.l? 


r 


V. In other cast's, quite possibly we may still have r \ . 
though the individual items mav be different. (Jenerallv in marriages 
X > // : in accidents 51 > : in vaccinations 51 < w : in total- 

abstainers N > // and 51 < m : in military enlistments and in many 
diseases N < a and 5F > m ; while in certain alarming epidemics, 
especially cholera, plague, and malaria, we find in addition that 
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I < i and E > e —so that in these cases, which are particularly 

*■ * 

our present subject of study, we generally have r > E. In fatal 
accidents M = 1 and X. I. and E all = 0. When considering 
happenings among the same individuals, we also put X. I. and Eat 
zero, though may be anything : and in otlier special cases we 
omit various elements. If the surrounding population is not 
affected, I — 0: and if the affected individuals cannot move*, 
E b • so that it will be seen that the equations can lie made to 
cover a wide fie ld. The theory will sometimes also apply to inani¬ 
mate objects as. for instance, m commerce, where the variation- 
elements may be taken to mean manufacture, waste, exportation. 

and importation, though the original equations may have to be 
somewhat modified. 

(iii) The i eve)sion-eleineut^ may also at present lie taken as a constant 
and be calculated for the same unit of time as is used for the other 
elements. In the case of independent happenings (to he defined 
presently) rdt means merely the proportion of affected individuals 

which become, in element of time, capable of being re-affected_as 

b} di'voice in maiiiage, and by discharge in certain employments. 
In dependent happenings, however, it implies also the loss of capacity 
foi affecting others. Thus, in infectious diseases (which are dependent 
happenings), it implies loss both of immunity and of infectivity— 
not recovery from sickness, which is merely an episode of affectedness 
fiom disease. In some diseases, such probably as leprosy and organic 
diseases, r is zero or nearly zero ; in others with long-continued 
immunity, such as many zymotic diseases, it is low ; and in others 
again, with comparatively quick- loss of immunity and infectivity 
such possibly as uasal catarrh or de„ s „e, it must be large. In 
yany diseases, however, it is quite unknown at present ancl must in 
act be calculated fronTthe integnatedequations (which, it is hoped, 
will prove of use for this very purpose). In the cisTbf stight accidents 
it is unity ; in fatal accidents it is zero ; and in snake-bite or heat¬ 
stroke where recovery is quick if it occurs, it is high, though at the 
same ime t le death-rate M is also high. In controversial parties 
due to rational divergence of opinion it should be high ; in ordinary 
party-politics, it is, in fact, very low. 
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(iv) Tl„* most important element is the happening-element, h. 
\\ e should clearly understand that in most cases the happening, 
whatever it is, does not select only the non-affected. but tends to 
lull on both groups alike. If, however, it chances to fall upon 
individuals who are already affected, it merelv re-affects them and does 

*/ * *J 

not cause them to pass from one group to another. Really the total 
number of individuals to whom the happening occurs in element 
ol time is hit . I\ that is, hit . (A A Z). But the individuals 
numbering hit . Z do not count because they are already affected, 
and therefore do not appear in the equations. On the other hand, the 
actual number of netr cases (which we may denote by F (It) is hit . A : 
a-Jid this quantity does appear in the equations 2 and 3. Ol* we may 
obtain F independently from the proportion 


F dt 


hit , 1 


A 

P 


( 5 ) 


Z). This is an important sub- 


so that, again. F =- hA — h (P — 
function in all cases. 

(v) Different kinds of happenings may be separated into two 
classes, namely, (a) those in which the frequency of the happening 
is independent of the number of individuals already affected ; and 
(b) those in which the frequency of the happening depends upon this 
quantity. To class (a) belong such happenings as many kinds ol 
accidents and non-infections diseases due to causes which operate, 
so to speak, from outside : and to class ( b ) belong infectious diseases, 
membership of societies and sects with propagandas, trade-unions, 
political parties, etc., due to propagation from within, that is, from 
individual to individual. In the former case, h or F will be constants ; 
in the latter case they will be functions of Z. 

-we inav 


To proceed now to the integration of equations J, 2 and 3- 
observe that the solutions, though tliev belong to the same class of 

* r * f 

functions, differ specifically in different cases. J he integrals are 
easy to obtain by ordinary methods and are familiar in connection 

ml V ^ -- _ 

with many statistical, chemical and physical problems : but the most 
elegant and useful method is to put all the cases into similar forms 
which can be brought bv substitutions into the same immediately - 
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integrable differential equation of the .simplest tvpe. We will first 
take a rapid survey of each solution as it is obtained. 


iSliUTIOX I\ 



penmgs. —Here/i or Pis supposed to bo constant. 
If h is^constant, the happening falls on the same proportion, hdt, 
of the population in every element of time. 

Z for A, using r and V for the variation-elements 


Putting P 


as in 4, and setting x = Z/P, equations 1 and 3 become 


- w ( W-i) ■* * 

dPjdt — eP — (t 


V) rP. 


X 


, ^ , *Tk< C-M P) -KV-<M~-r')nP 

dxP/dt = IP (1 — x) -f (V 

_ _ _^ ’ P "f* 

and, of course, = k.p ^ ) 


r 


/ 


) .rP, 


(<>) 

(7) 




h 


i 


dxP/dt = xdPjdt + P faUlt. 


Eliminating dxPJdt and dPjdt from these, we find that P cancels 
out also, so that we have s fCi -,c/ +(V-n-<)„(> - K ^jt - / -+ 


dx/dt = 7? — (/* + v — V + N + r) x -j- (e — V) x-. 
This has one integral form if v — V. and another if y. 


(8) 


(ii) Ihe Equivariant Case. —Here the happening is such that it 
does not affect the sum of the variation elements of the affected 
group though, as already mentioned, the separate variation-ele¬ 
ments may be different; that is, the total population is not altered. 
An example is, conversion to some philosophical creed without 
propagandists Another example is. the happening of slight accidents 
T" m wlllcI b however, r is generally unity;. and a third example is 
the happening of a certain standard of wealth, which, we assume 
tends to diminish the birth-rate, death-rate and immigration and 
emigration of the affected, each by an equal decrement, simul- 

... j In this last case, however, as the children of the affected 
will tend to be born affected, we should remove N from equation 2 
to equation 3 and omit it from equations 7 and 8. In all these cases 

v = V. 
























14 


Ronald Ross 

v - Y 

(iii) To integrate equation 8, however, write it in the form 


dxjdt 


K(L 


where Iv — h N d - 1 } md L 
il.rldf. vanishes when x — 


a. n 


Put y EE L — x ; so that x = L 


V 





dx I dt 


— dy/dt 
1 jy . dy 



(9) 


h'd v- 



Kd/. 


( 10 ) 


Intear atin a both sides of this, 


logy 


Ki -T constant 


If w 0 is the value of y at the beginning of the happening, 


log y 0 = constant, 


and 


y = y» e 


— K t 


•*<& 


( 11 ) 


Tlrerefore 


(L — x) = (L — x 0 )e 


— Kf 


and 


x = L — (L — x 0 )s 


— Kf 


( 12 ) 


This gives the proportion, x, of the total population A who are 

affected at the time t, this proportion being x 0 when t = 0. 

To find P. put v — V = 0 in equation 6, and we obtain a t onc e a 


i 


J U ILXL^l j. , ^ ~ J- 

inferential equation of the same form as equation 9 : so that 


P = Po**, 


(13) 


wliere P„ is the total population at the beginning of the happening, 
when t 0. This is an important function which expresses the 
natural increase of the whole population due only to the natural 
variation-elenKMitsT?.** "’Of"course v is small when the increase of popu¬ 
lation is not very rapid. When v = 0, the population remains 
constant, the births, deaths, immigrations and emigrations annulling 

each other. 


To find Z we have Z — xP. 

[n ifltiv), we defined F to be the number of ne\ 
that F = h (P — Zb‘ A If f = F IF , we have, when x 

>r ’~fr h (1 _ *) = h (1 _ L + Le _K< ). = l+L ' KC 

<-.KK>-V w ) 


ne'w~cases ambsliowed 
= 



( 14 ) 
































A priori Pathomefr/f I (IV) 


15 


1 lie functions for A and i/^ive what may be called respectively tie* 
(icfuxtl aiu\ proportional ntrres of affect ed indie/duals : and those for 
F and /five respectively the actual and proportional farces of w ,r 


cases. 


lo find / for any assigned value of ./■, we have 
K/ -■= log, (L — ./•„) — log, (L — ./■) 


2• 3025851 logj () [(L, — ./-„)/(L - ./■) |. (1,5) 




r 


The following equations are required for dnaKnim? the curves 

■ " '-■NtP -*• 

: «->- ,r .1 


and f : 


A (-y) = 

,= - & X - 

(&> c k C'H? L- f C. j - 

df/dt 


d-jjdt 1 



ty j.f l<o, ■JZJS'-I ~ k.(~k.y) 

- ?)• * (16) 
A-( /<• y ) 

x). ( 17 ) 

( 18 ) 

(iv) In order to analyse the cur^x (of affected individuals) we 
irst observe that as h, N, r are always positive, K and L are also 
always positive and L is always less than unity, while s~ KI diminishes 
and finally vanishes* as t increases. Hence' when, / is very lar-e 
■r inches the.limit L (equation-12) and never exceeds itV.r 0 ^°L<- 
Its Tangential (equation 9)' begins at the value KL (= h) when 
/o-^d^s remains positive, and gradually diminishes to zero 

" m [>~second tangential (equation 16) is always 
negative Hence as / increases from zero to infinity, x also increases, *' 

a rv,„Mnh1 increme',*, until it .pprotd,nates to the 

‘ w . 5 it is always conv^c'to the axis of t 

Ort^Wtan, curve / (e^JSon*^. of w o aS e.s) l,e., ills 

K,e “ feft value 4 = 0) and constantly di.ninishes'its 

IX ‘ 17 > '«“* •l'™™ negative. Its secon, 1 tangen- 


tial is always positive, and th^cufe^is, therefore, cbnvox to “the 
that is i\“ d ',', hen ,'. lS lar « e ^pproaches its final value h(l — L) = l 
be denoted by U^u\ ? ^ Let thlS COnstant 

happenings the ' , tll0U ®^ tbat; at Ieast 111 tlie ca « e of constant 

has been shown that this is not *, . ’ Dut lt; 

lot the case, because x cannot exceed 

































10 


Ronald Ross 


~£t t 


L. which is a proper fraction. What is it then that limits the value 


of ,r { 


Equation 8 may be written (when v = V) 
dx/dt = Ji (1 — x) — (N -f- r) x =f — (N + r) x 


(N -Vj) x. When it vanished 


(19) 



x = L 

ion of affected individuals comes to a 


and dx/dt vanishes when / 
and/ — l ; that is, the pro 

limit when the proportion of -hew cases exactly balances the propor¬ 
tion of recoveries-dind birth^s (it being assumed by supposition that 
the progeny of the affected are born not affected). If, however, 
X and /• are zero (and only then), L = 1, and the whole population 
does become affected in time. It should also be observed thaf® 



L) s ’JA./P , where A L is the final niniibei^oFnon-affected 
individuals when x = L and t is very large ; that is to say, even 
after tire limit or x has been nearly approached, the happening 
' continues and constantly affects the same proportion of the non- 
affected population as before—but now A also comes to a limit, 
namely P (1^~JL). 

! The easiest waT to assess the increase of x as a function of t is to 

* K 

o*ive to s k/ the successive values 1, 2, 3, 4 ... 10, 100, 1000 .... Then 

= 0, x will have the corre- 


it is clear from equation 



at, if x 0 


spending values 0, vL, §L, £-L, ... 0*9L, 0-99L, 0*999L, ... and t 

J -—-^ _ __u 


will have the values 


2-302 ... 

K 


log 10 (1, 2, 3, ... 10 , 100, 1000 ...). 


■Ujl =■ 


4n imjHtrtant point is reached when s 


K/ 




= 10, for then 

2-3 


f 


l 


K 


2*302 ... — t (say). 




At this point x lias increased to 9/lOths of its ultimate value L that 
is, most of the change of which is it capable has been effected, and it 
increases onlv slightly afterwards, namelv only by 0*09 and 0*00!C 

l » * ■ 

when l and 3t, and so on. On the other hand, a ieaclu s 


lato half its ultimate value very quickly when / 


i - 

A ' 


lie 


arlv. 


v 


% 

However large / may be, s _K ' never quite vanishes, and therefore 


x never quite reaches L. It is convenient, then, to find the value of 


. t H 
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t when, the number of affected individuals reaches the limit less one 
individual —that is, when Z = LP — 1. This is obtained at once 
from equation 15, since L/(L — x) = LP/(LP — Z) ; and the 
required value is 


t 


- log 10 LP 


(say) 


( 21 ) 


For example, if L = and P = 2,000,000 and remains nearly con¬ 
stant, the number of affected individuals is 900,000 when t — 
and is 999,999 when t = 6t. 

(The ratio xjt is the tangent of the radius-vector from the origin 
It. - the point (x, t)J If x () = 0 the value of this tangent is h, that is, 

dx^/dt, when t = 0 ; and is 9L/10t, that is, 0 • 390867? yvfien / = T ; 
- that, as the curve of x increases from zero to niire-tenths of its 

ultimate value, it Ces wfiolly vitlim the angle formed by these vectors. 

And the magnitude of tins angle depends only upon ji, and not upon 
the other elements, N and r. 

The constants L, l, and r may be written in detail, 

^ J 


L 


(Af+vf L ~ 


1 -f- (N" -j- /’) h 1 


Jr 1 + (N + r) 


-i * 


2*802 ... 
h + N r * 


( 22 ) 


From these equations it will bh easily seen that if h increases while 
N and r remain constant, L andu will increase and t will diminish ; 

but if h remains constant while\N or r or both increase, L and t 
will diminish and l mil increas 



1 and N ^ 0, l > L, since l 


■■—just as may be expected. If 

(N + r) L ; that is, as all the cases 
recover as soon as they occur, the total number of cases will not 

exceed the new cases during unit of time—as happens with slight 

accidents such as mosquito bites, where there is nothing in one 

happening to prevent a second or third occurrence of it If N 4- r = 0 

L = ’ - ’ - ^ - 5 


1, l — 0, and r 

1-151... 


h, and if h = N + 


1 7 

2? * 


hh, and, 


(v) It remains to examine some questions which require for their 
solution the integrals of P, Z, A, and F— all easily obtained. 


c 


- g-e. 
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Ironi equation 13 we have at once 


P (It 


P 0 £ rt dt = Vjv + constant 


rt 


0 


p dt = (P - P„)/y. 


(23) 


To lind the integrals of Z, A, and F, we integrate both sides of 


equation 7 as there written, or as written in the form 


and obtain 


dZ/dt = hV - (K - V) Z, {a 2 --j(tt 


to 




~ S -t c 


(Z'V) 


Z 


or 


F dt + (V 


N — r) Z cfe, «t z* 


(24) 


whence 




. 


A |P& —(K-V) Zdt; + 2 e 

h ^ - <K~V) " -< 2. 

P-Po 


(25) 


Z rf* = J 




v 


t Z — Z 0 
J A 


126 ) 


'/ p p 7 7 

Adi = (1 - J)h^A_L« _|_ Z- L 

* 0 ^ V 


0 


h 


(27) 


where 


J — kf( K — V) and'. 1 — J 


V —N 


r 


K — V 


And. of course. 


rt 


0 


rt 


F dt = h A dt 


* 


The same equations can be readily obtained by direct integration 

~ . ■£.. i *. e *P 

from equations 12, 13, etc. 

The concrete interpretation of these definite integrals of P, Z, and 
A, is that they express the total number of time-units lived during 
the period t by the whole population, by the affected population, and 
by the non-affected population respectively. 

For example, suppose that a population P 0 has doubled itself in 
1000 time-units, then, since lOOOu = log, 2, 4 v = 0-00069315, and 
the total number of time-units lived by all the individuals together 
o (2 — 1 )jv —that is, 1442*7 P 0 time-units. 

It should be observed that if v is very small, each side of 
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equation 23 takes an indeterminate form, the value of which is i ne 
value of P 0 (e v< — l)jv when v = 0 -namely, ;P 0 . That is, when tin- 
population remains nearly constant, the total number of time-units 
lived during t is tP 0 , as we should expect. 

The same integrals divided by t will give the wean numbers ol 

the total, the affected, and the non-affected populations respectively 
during the period t. Thus, if the happening is such that it does not 
affect the normal rate of increase or decrease of population, then 

the mean population will be (P — P 0 )/(l(5§ e P/P 0 )—as we should 
expect. <q. 

When multiplied by the appropriate constants, these integrals 
will also give the total number of the variation-events, births, deaths, 
immigrations and emigrations. Thus the total progeny of the non- 


affected will number n 


o 


A dt, and the total deaths among the affected 


t 

will number M Z dt The sum of the natural variation-events in 


a 


the whole population is (n — m + i 


t 


e) P dt , that is, P 

-o 


Po- 


Ldt gives the total reversions : and the integrals 

jr 


Similarly r 

multiplied by h give the total number of time-units lived in which 
the happenings have occurred—which is the same as the number of 
individuals to which the happenings have occurred. Thus 

't 

*J 0 A# 1S the total number of happenings (not necessarily for the 

first time) among the non-affected—that is the total number of 

new cases, P dt. —f.f*. 

* 0 

\It is now easy to interpret equations 24, 25, 26, 27. Since 

V — N — r = I — E — M — r, equation 24 means that Z, the total 
number of affected individuals remaining alive at the time t, is equal 
t e sum of all the new cases plus the affected immi gratio ns and 
^ S effected emigrations, deaths, and recoveries. Also as 

~ ^ “h 1 + M -(- E I, equation 25 means that the same 
quantity Z is equal to the total happenings plus the affected im- ? 

ngrations, less the affected emigrations, deaths, and recoveries, 

C 2 
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and also h*s> the happenings among the affected—which, as 
>tamd m III (iv). do not count. 


SECTION V 




u) Ijtt ittlcit! Happenni(j$ : the General Case : — The 

nwth-'d <>t working employed for the equi variant case has been given 
at smn«> length 1 'irnusc it will serve also for the other cases. We 
ii"w j>r• ■ c <■ i‘i 1 tn integrate the general equation 8 , namely, 

'/ r ih — h — (h -j- o — V + N -j- r)x -j- Lv — Y)x 2 . 


Writ.- tlii s m the form 


dr jdt — K (L — x) (U — x ), 


(28) 


when* K -- r — V and Land L' are the two roots of the quadratic m x. 
Let J j - - x -- [i and L' = a + (L where a — (h -f- N + r + K)/2K 
and 


fj 2 _ _ 


y? — hj K. These roots are always real because 


(h 


X -g / 


K ) 2 - 4/dv 

- (h - K ) 2 + 2 (h + K)(N + r) + (N + r) 2 , 


of which the left side is positive when K is negative, and the right 
side is positive when K is positive. And the roots are both positive 
when n > V. Observe that dxjdt vanishes when x = L or L', if 
this over occurs. 

0) Put 5 gg*~ z 


Then 


y = a/ ■ 

r/x L' 




(L' - my - 1 ). (29) 


vt-l tfr * 


ic 




1 . . I 


(!/ “ l ) 2 * 

hjjdt — K (1 


L dy /t 
• = Jvy (L 



Kyf 


L - L'V 


2 foe 


V 2 / - 

K (1/ - L) < 


(30) 


as in equation 11 . Therefore, as L — L — 2(3, 


L' 


X 


U 


x 




L 


L 


0 

” o 


X 


; ... 1 ) 


0 


i.-* ;«■ 


L-“tt 


or 




L 


2(3 (L — a: 0 ) 


LI 


(L' - x 0 y K » - (L - x 0 ) 

L' - L 4, '^f ^ 

- when 0 


0 . 


(32) 


/ -7 3 % 7 
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The first tangential of this is given in equation 28 ; the second 
tangential is 


- H — k — KC w < 2-n ^ £< v (_ L — nJ 

'dP = — 2K 2 (a — a) (L - a) (1/ • 


')• 


(33) 


(iii) The analysis of x, the proportional curve of affected individuals , 
depends upon whether K, that is y — V , is positive or negative. 
If it is positive, we are generally concerned with cases of Ini ui'ious 
Happenings ; because, as explained in Section III (ii), the nativity, 
mortality, immigration and emigration are likely to be respectively 
less, greater, less, and greater among the affected in such cases than 
among the non-affected. Non-infective diseases should be examples 
in point. On the other hand, if v — V is negative, we shall be gener¬ 
ally concerned with cases of Benefi cial Happenings, which improve 
the natural variation elements of the affected—such for example as 
conversions to total abstinence are said to be (if they do not diminish 

i and increase e in places where they prevail, such as the Prohibition 
States in America). 

If K is positive, e 2K # in equations 31 and 32 increases indefinitely 
with t (since p is the positive value of the radicle in the roots L and 
L ) , so that x ap proaches the limit L, that is, a — J3 * and never 
exceeds it, and therefore never reaches the greater root I/, that is, 
a + P* Therefore dxjdt is always positive for values of x between 
zero and L , and d 2 x/dt 2 is always negative for the same values, since 

“ p. Hence in this case x has a form similar to that which 
it ha^in the equivariant case. That is, it begins (when x 0 = 0) at 
zero aPap angle of which the tangent is h ; constantly increases with 
a decreasin^iqorement, so that it always remains concave to the axis 
of t , and finally^approaches the limit L when t is very large. 

When K_is nogative^tfiejquadratic expression in equations 8 and 
28 has only one change oP^i^n, and therefore only one positive 
root, L or a + (3. The pthS? ‘rbofi L, does not concern us because 
we consider only positive values of\ln this case, c 2K * diminishes 
indefinitely as t increases, and x finally^rhaphes the limit I/. Here, 
too, dxjdt always remains positive for possM&^alues of x 9 that is, 
while t varies from zero to positive infinity ; so that^tyagain always 
increases from zero to L\ If h + v - Y + N + r, thaHla^ 2Ka, 



a 
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<*- <- 0 


is positive. (Px/dt 2 is always negative in this case also, and x is there- 

i 

lure concave to the axis of x as before. But if this expression is 
negative. d*.r l df‘ 1 will be positive at first, but will change sign when 
;c —a. and finally vanish when x = 1/ ; so that, in this case, x will 
be first convex and then concave to the axis of t. This is the only 

k/ 

difference in general form. It is seldom that V — v can be greater 
than h X -j- r > which is always positive ; hence x has the same 
form in most cases where h is a constant. 


(iv) For the proportional curve of new cases,/, we have 


- x ); 


dfjdt 


1) dx/dt, 



h d 2 x Idl 2 . 


(34) 

(35) 


Hence, when r„ = 0, j begins at the value h and. as t increases, 
gradually falls to the limit h (1 - L) when lv is positive, and to the 
limit h (1 — L') when K is negative. The curvature is the opposite 

of that of x —as in the eqnivariant case. 

To find the value of P. the total population, write equation 6 in 

the form. 


Therefore 


1 /P . cZP /dt — v — Iva:, 


log P 


vt 


Kx dt 


vt — KLi 


2K(3s 




IK fit 


~ ~-K 


= (v — KL) t + log (y 0 — £ _ ' Kw ) + constant. 

On evaluating the constant we find that we can arrange P in either of 
two ways, namely, 


P = PoS 


vt t/o 


„-2 Kpt 

- e~ KU or P 


Vo 


1 


p -r< ?/o s 

- 1 - 0 W j 

Vo - 1 


>Ket i 

r _ — KL7 /•: 


. (36) 


If x 0 — 0 and t is large, these become, when Iv > 0 and K < 0, 
respectively, 


P 


p -Vt 


L' 


0 


L'-L 


s -KT -' or P 


P- 


L 


rt --K1 :t 

A f ^ ■ I 

0 I/-L 


(37) 


The former expressions are suitable when Iv is positive, and the 
latter when it is negative. Thus in the case of injurious happenings 
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when K is positive, the population will diminish indefinitely unless 
the natural increase denoted by v is large enough to compensate for 
the decrease, KL. due to the happening (L being positive in this case). 
In the case of beneficial happenings when K is negative, L is also 
negative while L' remains positive, so that the happening enhances 
the natural increase of the population due to the happening. 

As with equivariant happenings, x never quite reaches either L 
or L' ; and it will therefore be useful to find the value of t when the 
affected population equals the limit less one individual, that is, 
when Z = LP — I or I/P — 1, or when x = L — 1 /P, or 1/ — 1 /P. 
This can easily be ascertained from equations 29 and 30, and we 
have to find t from the equation 1/P = (L' — L) j (y — 1) when 
K is positive, and from the equation l/P = y (1/ — L) / (y —1) when 
K is negative. From these we have when x Q = 0 


log, LP 0 


(KL' -v)t = 

(KL — v) t = log £ L'P 0 


(K > 0), 

(K < 0). 


(38) 


To find the integral of P, we write P in the form (when x 


o 


0), 


P 


L'P 


0 


L' -L 


(y-KL) t 


LP 0 


Therefore 


L' — L 


J«-KL')f 


0 


P dt = P 0 e 


zt 


L' 


-KU 


L' -L 


v 


KL 


t T 

& _ P r.vt 

- r 0 z T1 


— K L’t _ ^ — vt 


L'-L 


v 


kit • (39 > 


The integrals of Z and F are found from this by the same methods 
as were used in equations 24 and 25, or by integrating both sides of 
equation 6. From the latter we have 



V) 


■t 


0 


Z dt 


v 


t 


0 


P dt 


(P - Po). 


The remarks made in Section IV (v) apply here also. 



Section VI 

Independent Happenings .■ Constant New Cases. —A third example 
of independent happenings remains to be considered, namely, the 
one referred to at the beginning of Section IV, in which F, the number 
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of new cases, is always a constant. The only instances of this which 
I can recall are the cases of certain societies (such as the Royal Society) 
T\hich elect a fixed number of new members every year. Indeed, 
this process of increment seems to be too artificial to be seen in 
nature, unless, perhaps, in certain cases of seeding or spore formation, 

and we need only note the form of the principal function for future 
reference if required. 

Here, since F is a constant, let F = 
become 


c. Then equations 6 and 7 


dP/dt = vP - (v - V) Z, (41) 

dZfdt — c + (V — N — r) Z, (42) 

= K (L - Z). 


where K = M — T -)— E -j— r and L = c/K. The solution of this is 
given in equation 12, except that Z is here substituted for x. If 
Z 0 = 0, we have 



1 g— (M — I+E— 

C M-I + E+r' 


It is to be hoped that the case of the Royal Society is at least an 
equivariant one. 


Section VII 


(i) Dependent Happenings: Proportional Happening .—Referring 
to the definition of dependent happenings in Section III (v), we see 
that in such Ti must be a function of Z, and consequently of t. First 
consider the case in which each affected individual affects or infects 
c other individuals in unit of time, c being a constant. This may be 
taken as being a first approximation to the study of those very 
important happenings, the infectious diseases. 

The total number of the happenings which occur in element of 
time will be cZdt ; but, as in the preceding cases, some of these 
may chance upon individuals who are already affected, and the 
number of new cases in element of time, namely, Fdt , will therefore 
most probably be given by the proportion 


FdtjcZdt = A/P, 



that is F = cZ (1 — x ), and h — cx. 
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(ii) Equations 1 and 3 now become 


cTP/dt = vP - (v - V) xP, 
dxPjdt = cx P (1 — x) —(- (V — N — r) xP 


(44) 

(45) 


Treating these as in Section IV (i), we find that, as before, on 
eliminating dxPJdt and dPjdt, P also cancels out, and we have 



(c — v + V — N 
Kx (L — x), 


r) x — (c — v -f V) x 2 , (46) 

(47) 


where K = c — v -p V and L — 1 — (N -j- r)j K. Here again dxjdt 
vanishes when x — L. Integrating by substitution as before, put 


Hence 


y = (L — x)jx 


x = L/(l + y). 


(48) 


L dy 


(1 + yf dt 

dyly 


Iv 


L 


L 


x 


1 + y \ 

i + y/ 

-KL dt, 

y — y»* 

L 


1 | (L J'Xq 

1) s _KU 

Lx 0 z KLt 


a? 0 s KU + (L - 

-*o)‘ 


IvL 2 


y 


{\-Pyf 1 


- K U 


(49) 


(50) 


Befoie analysing this it is advisable to obtain the integrated 

expressions for P and Z. Dividing equation (44) by P and integrating 
both sides, we have 


log P 


vt 


(v-Y) 


xdt 


whence 


and 


vt 


P 


(v Y)/K . log (x 0 e KLt -{- L — x 0 ) -f- constant, 


P - rt 


L 


(c-V)/K 


0 




KLt 


+ L 


x 


0 


(61) 


Z = Z 0 e rt £ KI ' 


L 


Cl K 


XqE 


KLt 


+ L 


X 


o 


( 52 ) 


where the fraction within the brackets obviously = (x/x 0 ) 


-KU 
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I !:>• • • x ! > i i< • ir \':i| iii' 11 1 / |< ..jvcn by 


KlV - loir. --^ 


.r 


loir. 


L 


x 


0 


X 


(53) 


tiiii >r analysing the proportioned cnr> s e of affected individuals , x , 
tii-. iit'i r.tnx-Mtial is iriven in equation 47. The second and third 

t a 1 1 g> • 11 11 . 1 1 - ai’i* 


il'.i dt~ K-.r (L — (L 


x) 


dr-' 


K 3 ./ 


( L — x) (Lr — 6Lx + 0.x 2 ) 


(54) 

(55) 


J' < »I I 1 t 

* 

_»' 11 i 1;: 

c m r \ * 


i ;1 1 1 * [ 1 30j it Will be SG 611 that till tll6 tail- 

<>r .r jin- /.fin when .i -- u and = L. If x 0 is very small, the 

..-•-ins at nearly zero when / (); and then increases very 

* 


slowiv at tir.'t and more rapidly afterwards, remaining convex to 
it,., axis of / to begin with. When, however, x = -|-L, there is a 
eliaiii;i. nf curvature and .r becomes concave to the axis of t, and then 
tinally approaches tie* ultimate limit L when t is very large, and never 
.,xca*ds that limit. Its maximum rate of increase, namely ^KL 2 , 
is reaclu'd wlnm ./■ — ITj ; attd then 

I sAj — log [(L — x 0 )/x 0 ]. 

1 lie curve is now seen to be a symmetrical one with this point as the 
centre of symmetry. Moving the origin to the point by putting 
s ----- \L + x' and t = 1/KL . log [(L — x„)/x 0 ] + t\ we obtain 


x 


ux 


„-IvT X\ l/i i - -KI/ 


)/(! + 


)• 


(56) 


Substituting — l* for t ( in this, wo find that x also merely changes 
its sign wit!iout changing its numerical value. In fact, x has the 
general shape of a long-drawn-out letter S. 

-I-. -■ r rt 

(iv) The curve dxjdt is also a symmetrical one with its centre ol 
symmetry at the same distance from the primary origin, namely 
Jog \{L — x„)/x 0 ]. When x 0 is small, the curve increases slowly at 
first, then rapidly; reaches a maximum value, |-KL 2 , at the centre 
of symmetry, and then falls just as it rose. It has two changes of 
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curvature which are the two roots of the quadratic 
d 3 xjdt 3 , namely when ,r = (-1 \ ) L : that is when 


expression 



x' = T0-28867L, 


and IvLi' 


T log [(\/3 



- 1)] = =F 1-31697. 



At these points the ordinate of dxjdt , that is, of Kr (L — x), is 

KL 2 + \Z-fy) (1 — | — VlV)> 

that is § jKL 2 , or two-thirds of the maximum ordinate at the centre 
of symmetry ; and the ordinate is, of course, the same at both points 

of change of its curvature, as also at other points equidistant from 
the centre of symmetry. 

The curve of dxjdt is therefore a Tegular bell-shaped curve, very 
similar to those often found in epidemics. 

By the conditions of this kind of happening, x 0 and Z 0 can never 
be zero, for if they were there would be no new cases. It is con¬ 
venient therefore to take Z 0 as being one individual—that is, 
Z 0 = 1, x 0 = 1 /P 0 , and dx 0 jdt = Iv (LP 0 — 1) /P 0 2 , or nearly KL/P 0 . 

As with the kinds of happening previously considered (see 
equation 21), x and Z never quite reach the limits L and LP re¬ 
spectively ; it is useful then to find when the number of affected 
individuals reaches the limit less one individual, that is when 
u = BP 1 and x = L — 1/P- If Z 0 = 1 and x 0 — 1/P 0 as just 
suggested, then, owing to the symmetry of the curves, this figure 
will be reached when x = L — x n ; that is (equation 53) when 

KLt = log [(L — x 0 )jx 0 ] — log 0„/(L — a; 0 )], 

= 2 l°g (BP 0 — 1) = 2KLi:(say). (58) 

Thus if the total population is large and L is not very small, t, the 
abscissa of the centre of symmetry when x 0 = 1 /P 0 , is very nearly 
equal to 1/KL . log LP 0 ; and at twice this period from the beginning 
the limit of affected individuals less one will be reached. We 
observe also that when t - 2r the ordinate of dxjdt will be the same 
as it was at the beginning, namely, nearly KL/P 0 . 
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If x 0 — 1/P 0 , the tangent of the radius-vector from the origin to 
the summit of dxjdt is |KL 2 /t, that is 

K 2 L 3 /4 log (LP 0 — 1 ). ( 59 ) 

From equations 48 and 49 we have 

-ku _ x_ L :f *n aP LP 0 a -pPp Z LP 0 — Z 0 

.r 0 L — x *oPo LP - xP ~ Z 0 LP - Z ‘ 1 ; 

this provides a useful formula for a cursory estimation of the 
increase of Z with respect to t when the latter is small. For if P 0 and 
P, Z 0 and Z are not very different, and LP 0 is large and Z 0 small— 

as at the commencement of the happening—-then this equation is 
nearly 

Z - Z 0 e KL{ . 


Suppose that Z 0 is one individual (t — 0), then Z = 2, 3, 4, ... 
individuals as s KL( increases through the values 2, 3, 4, ...; 
that is, as t increases through the values 1 /KL . log 10 * (2, 3, 4, ...). 
This same thing does not hold true, however, when Z is larger. 

Similarly, we may have at first F = cZ 0 e KU — e (Z 0 e KL *) 2 . 

to 

(v) The actual and proportional curves of new cases are respectively : 



F 

- C Z (1 

— x), 



f" 

= cx (1 

— x). 

(61) 

To analyse f 

we have 




dfjdt - 

— cKx (L — 

■ *) (1 - 

■ 2 *), 

(62) 

d z fjdt 2 = 

= cK 2 x(L - 

x){L- 

- 2 ( 2 L + l)x -F 6 a; 2 }. 

(63) 

And, since L = 

— 1 — (N -f r) /K, we may also write / in 

the form 


K//c 

~ dxjdt ~J” (N -j“ x. 

(64) 


The form of the curve / depends upon the value of L. If L = 1 , 
that is, if (N + r)/K is zero or very small, the curve / equals or 
approximates to the curve dxjdt multiplied by the constant c/K 


* As in original, but may be log* 
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which, is generally greater than unity. That is, in this case / is 
symmetrical bell-shaped curve with the properties just described. 
In any case, however,/is such a curve plus the constantly-increasing 
function c(N + r)xJK, or c(l — L)x. We may, therefore, expect 
that, at first, when x is small /will follow the graph of clxjdt, but that 
as x increases, the ordinates of /will become increasingly greater than 
those of dxjdt. At last, when t is very large and x approximates 
to L, dxjdt becomes nearly zero and / approaches the constant 
value cL(l — L), at which it remains indefinitely. Thus/is quite 
symmetrical only if L = 1 ; but is nearly symmetrical if L is only a 
little less than unity. 

It is evident from equation 62 that dfjdt vanishes when x equals 
either L or That is, if L > -J-, / reaches a maximum when x — V 
the ordinate of this maximum being exactly ci(l — or \c ; and 
as x increases above this value, / then decreases, and falls towards 
its ultimate constant value cL(l — L) as x approaches L. If 
L = i exactly, dfjdt vanishes only once (and that not quite), namely, 
when x = L = i so that / now always increases until it reaches its 
maximum, which is also its ultimate value. If L <g however 
x never attains the value and consequently dfjdt has no vanishing 
point for any values of x between x 0 and L, and therefore again 
always increases for all considered values of x and therefore of t 
Hence the curve of/has two forms : (Type I) if L > i an irregular 
bell-shaped form (becoming regular if L = I) ; and (Type II) if 

L > h a drawn-out S-shaped form somewhat similar to that of x 
itself, but not usually symmetrical. 

From equation 63 it will be seen that / may have two changes of 
curvature as a- varies from z 0 to L, namely, at the two roots of the 
quadratic expression in the value of d 2 xjdt~. The roots are 


6x 


(21, + 1) ± V[(2L + 1) 


6L] ; 


(65) 


and are real and positive. The lesser root is always less than L, and 
is therefore always attained by a; as it varies from z 0 to L—so’that 
/ always has at least one change of curvature. The greater root is 
also less than L when L > *; is equal to L when L = / but £ 
greater than L when L < so that the second change of curvature 
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occurs only when L > J. As x never quite reaches L, there is no 

second change of curvature when L = / These results are therefore 

as were to be expected from the two forms of the curve of / already 
discussed. 

As already seen, when L > / the summit of the/-curve is |c and 
is reached when x = The corresponding value of t is 

t = T — i/KL . log(2L — 1), (66) 

the logarithm being negative since 2 L — 1 , that is 1 — 2 (N + r)/K, 
is less than unity. 

It has also been seen that when L > / falls after reaching its 

maximum, \c. to a constant value, cL(l —L). The difference 
between the maximum and the ultimate values of / is therefore 

jG (2L — l) 2 , 

and their ratio is 

1/4L(1 — L) = l (say). (67) 

As 'with independent happenings (e.g., equation 22 ), we shall use 
the symbol l to denote the limit which / approaches when t is very 
large, and this will be examined further in Section VIII (iii) (Part II). 

(vi) The curves of new cases in proportional happenings are 
especially important because, if the original assumptions upon which 
they are based are sufficient in themselves to explain time-to-time 
variations in the frequency of infectious diseases, these curves should 
agree, at least for a first approximation to the truth, with the curves 
of new cases actually observed in such diseases. Now, according 
to Brownlee (Section I), the curves of many epidemics, especially of 
the short and sharp zymotic diseases, tend to be remarkably sym¬ 
metrical bell-shaped curves, roughly similar, in fact, to the /-curve 
developed from certain values of the constants. 

We now proceed, therefore, to examine these constants with greater 
care. They are in detail, 

K = c — v + V. KL = c — v + V — N — r, 

L = 1 — (N + r)/(c — v + V). (68) 
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'1 ho particular elements have boon alrearlv discussed in Sections 

4 ~- 

nr and V. \ ho natural-variation element v governs the natural 





increase of the population apart rrom exceptional nap 
but V, X. r are cor unacted with the .special happening under considera¬ 
tion. If the period covered by such happening- is short comj>ared 
with the average life of the individuals concerned, v may be so small 
as to be negligible, and X also may be verv small. Tn injurious 
happenings such as infectious diseases, v — V should lie positive 
'Section V (iii) ], and V will be negative if M > X and E > I. so 
*■ -' V still remains positive even when v = 0 ; but, as a rule, 
except in very' fatal diseases, o — V is very small. The element 

4s 

1 oiay be unity when reversion is immediate, as in slight accidents, 
but may be neai Ij zero when the minium tv conferreii b^ ' one attack 
ol a disease lasts lor a long time on the average, as in many zymotic 

\ and t are only 




V 


vs. 

r I 





SO 



^ ~r is usually a small or very small positive fraction in 
the case ol infectious diseases (its actual magnitude depending, of 
course, on the magnitude of the time-unit taken as well as on the 
elements themselves). The element of special importance is c, that 
is the average number of individuals infected or reinfected in unit 
of time bv each already infected individual. 

%r 




v, V. N, r 

may be taken as fixed for the kind of happening considered, c is 
generally unknown, and. indeed, one of the ultimate objects of such 

studies as these will be to ascertain its value. At present, however, 
we take it to be a constant. 

We observe fust from equation 53 that if KL is negative x must 
constantly diminish as t increases, which means that a change must 
ha\e occurred in at least one of the original constants. 



assume at present that KL is positive, and reserve the study of 

the case when it is negative for Section IX (Mvpometric Happening • 

Part II). 11 

Secondly, as KL is taken to be positive, c. whatever it is, must be 
greater than v — V 4- N 4- r —generally a small positive fraction : 


and therefore greater than each item individually. 
Thirdly, as c « - V + N -f r and X and r 


V 


\ , that is, Iv is always positive 


are positive. 
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Fourthly, as KL >0, K > N + r, and therefore L is always 
less than unity (as it must otherwise obviously be), and always 
positive. Clearly also, L increases with the increase of c and V, but 
diminishes with the increase of v, r, and N (the last of which exists 
in K but cancels out from KL). 

We must now examine the conditions which hold respectively if 
/is to be of Type II or Type I, or is to be nearly or quite symmetrical. 

It will be of Type II, that is, always increasing in a drawn- 
out S-shape [(v) above] if L is not greater than \; that is, if 
c } {) —. V ~p 2(N -f- r). But as just shown, in the cases which we 
are now considering, c must always be greater than v — V + N + r. 
Hence in order that / shall be of Type II, c must lie between the 
limits v-Y + 2(N + r) and v — V -|- (N + r). Now if N + r is 
small, as is usually the case in short and sharp epidemics of zymotic 
diseases with long immunity, these limits will be narrow ; in other 
words, out of the whole range of possible values of c, that is of in- 
fectivity, only a small sector, and that the lowest one, will make / 

of Type II. 

On the other hand, / will be of Type I, that is, of a bell-shape, 
regular or irregular, if L > ]—that is, if c is any number greater than 
v _ V + 2(N + r). In this case / rises to a maximum ordinate 

\c and then falls to a constant level cL(l — L). Lastly, / becomes 
nearly symmetrical (see (v) above) when L = 1 nearly, that is, when 
c __y__y_p(l -j- y) (N -)- r), where y is considerable ; and also 

of course when N and r are very small. These are probably just the 
conditions which hold in many of the short and sharp epidemics of 
zymotic diseases, such as measles, scarlatina, and dengue ; and they 
produce curves which differ from the perfectly symmetrical curve 
of c/K . dxjdt only by a small increasing term which, however, never 

exceeds c/K . (N + r) L. 

For an example, suppose the case of a human infectious sease in 
which « — Y = 0-002, N - 0-0004, r — 0-02, all assumed for a 
time-unit of one week. Then KL will not be positive unless 
c > 0-0224 ; that is, the happening will have no effect unless, say 
1000 affected individuals infect or re-infect 22 • 4 individuals per wee 
on the average. If 1000 such persons can infect or reinfect 42-8 
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persons per week or less, the /-curve will be of Tvpe II and will 
always ascend to the limit \r (that is. 0 0107 if r = </-042S), at which 
figure it will remain indefinitely. Tf. however, r is greater than this — 
if the 1000 affected persons can infect or re infect, sav. 100 persons per 
week then the /-curve will rise to a maximum value of 0-025 of 
the total population, and afterwards decline to the constant value of 
0-01047 of the total population, or 10-47 per thousand, at which it 
wdl remain (L being 0-792 in this case). If c is still greater, if each 
affected person can affect one other person per week, so that c 1, 
then L 0-9795, that is. nearly the whole population ultimately 

„.... rcaclira ; 

maximum of one quarter the whole population ami thereafter declines 
to the constant value of O-OIKOT, or nearly 2 per cent, of the popula¬ 
tion, at which ,t remains indefinitely. (For further examples set 
Section f\. Part II.) 

(vn) Integral expressions for P and Z were obtained in equations 

31 and o2 by integrating both sides of equation 44. These may be 
written J 


* 

j 


^ r # e‘‘ (x/x 9 . s Ku )' -'' K 2 = P 0 (xjx 0 y- V > 

Z M {x/xq y 1 k £ rt g - \ > u 


Z — Z 0 S H e klt (x J' 0 w-KCO«;K 


) 


(69) 

(70) 


Two important special cases present themselves here If y = (j 
the happening has no effect on the total number of the population,' 

’ "‘ l ‘ , " , “ ,10U3 d,seal « is *>®l, “» chicken pox or (?) den»ue or 

CL^Tu Thi “ ab0a,J ^ CMeH th ° ^uivariant Ca°se of 

Proportional Happenmg. In it, P - P„ e » as in Section 1V (u) and 

of coinw z = Z.,-. But if e _ v is email without being zero, care 
must be taken with the functions if c, N, and r are small. 

Caw L°twfTv CMe ,T' 1 ” Wien ” = °' ma y be “lied the 
2^ l mtant N t tU / **<**«■ It ie particularly important. 

.pp^.:r^ , ^ of ;h:n of : h : imi ;™ ,uaie c “ 

thTTmTv «th the other constants 

may be neglected; in fact, the a priori curves of most 

pidemio diseases may be considered on this Hama 
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If v = 0 and V is negative (as in all injurious happenings) we sec 


that, since 


X)-T3 / rp //y» \ — 1 

JL - JL q ^ 


(71) 


the population diminisl 
of course, what was 



s indefinitely as time progresses ; which is. 

since by supposition there is no 
natural increase to compensate for the loss due to the happening. 

Y = 0 then, since ./■ -- L ultimablv. 


If neither v = 0 nor 


v 


that is, a constant, the question whether the population will finally 

? sign of the expression 




increase or 




v 


( v— V) L = V -j- \{r — V) (N -|- r)\/(c — v -p V). (72) 


An important value of P is reached when the ./-curve attains its 
centre of symmetry, that is when x — JL and 


t 


1 /KL . log [(L •b>)/> / 'o]- 


(see (iii) above). Then 


P = P 0 (L/2.r ( ,r- v " K [(L - .r 0 )/r 0 ]' 


r — (r — V) I.I/KI, 


(73) 


P» (£) 





if x 0 is very small, as when it equals 1/P„. 

Another important value of P is reached when Z 
is when t 2x (equation 58). 


- LP - I. 




P = P« [(L — -r 0 )/.r 0 ] 


{2r — (P— V) L| K I- 


(74) 


A third important value of P is reached (if L > J) when the J -curve 
attains its maximum ; that is, when r = J and t — (equation Gib- 
Then 


P = P 0 (l/2.r„) 


(i--Vt.lv 


• r o), (-L - iy_ 


‘ [ r — I r — V ) f J hi- 


. (75) 


[ 1 ; . (Ij •* of, V — 

These expressions are required to estimate the actual values ot 
Z and F at these points ; hut. if L is not very small, t will not be very 
large [see Section VIII (vii)], and therefore P can vary but little in 
that period from its original P 0 . and, at least for rough calculations. 

can be taken as remaining constant. 

(viii) It remains to consider the integrals of P. Z and F. which will 
be required for the analysis of several questions. e proceed as in 
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.>'J 


fix- ca-i- of constant happenings |Section [\' (v) I hv integratin'.; both 

I | 1 

'!'l ,-> f -l"* fundamental differential equations (44 and 4b). Hence 


J 


P 


IJ 


/' 


u 


IV/ - (V - V) I /At 


(70) 


f J 


z z 


0 


Id// .. (v 


/ 


J IJ 


X — r) I /a 

v o 



(77) 


r r t . 

lo interpret the first of these equations, we note that the actual 
ehange in the total population effected during the period /. namelv. 

I l*„. equals what the change would have been In* natural- 

% 

\ dilution only, namely J el\//. |i*ss th<* dilhuvnei* between the 

. O 

natural-variation and the happening-variation in the affected part 

of the population, namely (/• V) '/At. The second equation. 

which is the same as equation 24. has been already interpreted in 
S*< tion l\ (\). ft means that the total change in the number of the 
affected .lurin', the period /. namely. Z-Z,„ is equal to the sum of new 

casts. namel\ | lo//. plus the affected immigrations, and less the 

* ( f 

affected emigrations, deaths, and recoveries, namelv. 


(f E - M 


(J 



Ifotli these interpretations agi(‘e with what was to l>e expected. 

T,1 “ int " yni1 ' " f F ^ he obtained at once in the important 
case when /' - o , and are 


- O 


/At - (1> - p n) 


( 



F.// - Z - Z 


0 


u 


(\ 


X _ r) (P - p„) \-. 


(79) 


If in these we have also \' _ u. then (P 

ate since p remains constant : but tile 
from equation 71. 


*• id ^ becomes indeternun- 
value of the fraction then is, 


< (J - K) \ 


Lt 


1 


Hie integrals of Z and F can also 

m 

integration ot their values when r = 


he easily obtained bv direct 


0. 
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We should note that the total births, deaths, immigrations, 
emulations, and recoveries among the affected during the period t 

aie gi\ cn, as with constant happenings, by multiplying VjI,1 hv 

i i 

* 0 

N, AL I, E, and r respectively. 

Tlie main result has already been summarised in the prefatorv 
section. We can see from (vi) above that, for manv values of the 
constants, the solution of the fundamental proposition stated in 
Section II yields curves generally similar to the curves fmpientlv 
found m epidemics that is, curves winch are bell-shaped and nearlv 
symmetrical, and tend to decline more slowly than thev rose. This 
decline is m the case of the ci prion curves due mere I v to the exlians- 

L 

tion of susceptible material ; but the further studies contained m 
the second part must be considered before attempts are made to 
apply the results to actual observations. 
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PARTS II and HI 


PREFACE 

Parti of this paper was published in the ''Proceedings,” 1916, 
A, vol. 92, having been read on November 11. 1915. In June 
? J 1 OClt t A was kind enough to give a Government 

Grant for providing me with assistance in order to complete the 
paper, and for carrying on further studies upon the subject ; and 
Miss Hilda P. Hudson, M.A., Sc. I)., was appointed for the work from 
May 1, 1916. The continuation of the paper has accordingly been 
written in conjunction with her ; and I should like to take the 

opportunity to express my obligations to her for her valuable assist¬ 
ance, especially in regard to Part III. 

The entire paper here presented is still limited to the theoretical 

side of the subject, as defined in the Introduction to Part I. 

Records of epidemics are now being examined in order to find how 

far the theoretical results which we have reached may be applied 

to them ; but these studies must be reserved entirely for future 
discussion. [They were never completedj 
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Section VITT 


(1) Hypothetical Epidemics: Deduction of Constants .—In any 
case of independent happening we can use conations 8 and 12 with 


V. e. X. r. all zero. 


Then 


d.ridt — hi 1 — ./•), 


and, if — 0 


ht — — lo<i f (i — 


x 


1 ...2 
kji ** ■ 


; !j,r 3 ... if x is small 


2*30 log 10 


1 


1 — ,/ 


for all values of x. 


If the value of x is known for any value of t . say a year, this equation 
enables us to find h, and tlience x for any other period, hour, day, week, 
etc. 

Thus suppose that child-birth, death, or migration occurs to 20 
per mille of a population of all ages during one year, to what propor¬ 
tion do thev occur in one day t W e have 

365 h = 0*02 + 0*0002 + 0*0000027 + . 

h — 0*000055. 

To find the proportion for one hour, we divide by 24 : and to find 
it for 7 oi* 30 days, we multiply by these numbers respectively. 

Or suppose that 85 per cent, of plague cases, taken at any stage 
of the illness, die within 30 days, and we require the daily death-rate 
of such mixed cases. Then we have 


30 h 


2 * 30 log U) (1 — 0*85). 


h — 0*003, 

giving a constant death-rate of 6*3 per cent, per (hem. If the cases 
are all taken from the commencement of the illness, the curves are 
quite different, and show the highest death-rate some days after the 
onset—but this does not concern us at present. 

If the annual birth-rate is 30 per mille , and the annual death-rate 
is 20 per mille, then the corresponding daily rates are 


h = 0*000083 


and 


m = 0*000055. 
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if the iimnigflti' tt and emigration rates counterbalance 
«-ach other we nuy take fur the daily variation element 


n 


m 


l >-iHHHrjH 


In rhi^ example increase or ’ In* |><j|mlatioa at the end of the 
V'' 1 *! WDBld amount to 1 |» rent, that ix, I* would vary from I*, 

at the beginning of the year to |*d F 0 at the end of it. But if we 
uae the value of 9 just given lor equation 13, namely, F = F„; 
(putting t ■ 306 and vt 0*0103) we obtain I* — |-oi<>3 F 0 at 

*d th» •* ir -which in a little greater than the figure jint 
K ,VVh - rhe raaeon for this i- of that in equation 13 w< 

aaauane that rh. pro p a y tion ia eontmuoue—that in, that during 

the year, the pr. •-••n\ of the original population will them->elve 

have progeny. But this is not aseune <1 in the value or I' in -r 

»-orametl. 1 1 ‘J v • 




We may ilxi calculate «- from equation 13 if we know the values 
of F, ami P at the beginning and end of t. Thus the population of 
Inghnrl ami Wales wea I0JUM 236 in lui 1, and 3»» o7'»,4irj m i ut i 
Then Mnee at * log. F log. F, and I = 365 > I0o, tin daih 

v * r ** t * ul1 *• • — O tsgl0i47, the yearly variation element 

ia 0 01 -Jtfrtd, and the ten-yearly variation-element 0* 12000 If a 
Imputation doubles itself m a period t, we have « ' * log. j o « 93 , 
-o that at the above rate th Kuglish population ought to 

•5.5 years -which ia about the 

1 i" i ■ * i’t m I 1 a I | fat . s» nii,_ 



t'» < i y death-rate of plague is taken at the figure just calcu- 
iMol. of 0-OM, and the plague birth-rate is a normal one of 30 per 

mUU • J*' <*o that N a), the* V = N - M = - o«n;3 - . 0 

that., d r * O-UMM* r - V - 0 0C3 ... -lady. Such figures’give 

*** eomaete ideas of H fiat the C<»iat<uitit are likely to lie in cases 

of human ej*b-fi«e- (Compare Sections Iff and V (iii).'| 

The Ifc version-element r is particularly elasive since we seldom 

to ^ *** m^mAxnu It rt is beet to proceed 
“ W * >lMi “ *f mn * etion equation ^.4. and to assume tliat 90 

p*r cej af th- attoied individual, r-vert i„ the time /, *h-r 

Hfj » Mkor mure vears. . S b 
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1«* remind* d that KL m the coefficient of l ill -r, the 





nroDortion of affected individual-) (collation .30), ami that die /-curve 

if ' * ' 

i L, and its maximum value 


reaches it* centre of nymiiiefrv when / 

* m 


1, when t it v«f large lS*-«fion V II (iii)l. 

/ « i ' /1 


reach- - lift centr*- of avmmetrv 


The period at winch x 
and dr dt reaches its maximum is 











i 
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KL 


loo 


I. 


a 


* # 

./ ,1 


• t! f u - l/Hii and LP„ i- large compared 





with unity then as stl^ **ted m connection with equation .38, we 

may now w rit* 





i. I*, +■ !•% Y — <T r 1 ). 


tMiea thu |»-r i>m| tIn- maximum numlier of ,ifT. < f..I imhvuluals 


e^ oiie Is fear he* I, 







UoUrlv, when I it very Large, tie* ultimar* due of tie* proportion 




of newr ctflM lift/. But if 1- - \ / n->* ^ to a previous maximum value 

ft* 

ftWft f ■!, ami th*-r< after falln t<» ir <* ultimate valm- / I Section 
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rhe ratio bg|« n tin- maximum ami t!**- ultimate 
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* fj**rio<l at which th** maximum value 
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anti when ' 



r., 2L - I 



imt I.P. in large, this may he 
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But ttm maximum of/(pr«v ,.m, to it* ultirnat.- value) occurs only 

if Y > * 3 
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J'lien for a daily element 365 


r 


1 , 

t 


^•30. Hence if 90 per cent, of 


tlK ti ‘ VCar ’ the daily aversion-element is 

' ‘ ( ()3L 1Ilus tlle aversion-element is likely to be lame 

com pa led with r or X and correspondingly important, and will be 

n no nuan.s negligible even when the considered disease is one in 
winch the acquired immunity is supposed to be very lasting [com¬ 
pare (Section rrr (in) i. For a case reverts not only when it becomes 
again capable of showing a recognisable infection, but when it is 
again able to harbour the infective agents sufficiently to afford them 
a an/as whence they may infect others ; and it is quite possible and 
indeed probable that this may often occur much earlier than we 
iina.im . Ihus poisons vaccinated against smallpox may acquire 
t le disease m a very mild and modified form after only a few years, 
and may then infect lion-protected persons in full force ; and there is 
n j pi oof that a mild hist attack of measles or other diseases may not 
gi ve quite as short-lived an immunity—that is, that those who have had 
the disease in infancy may not comparatively early acquire it again in 
listinguishable foim and then spread it. In such cases, r maybe 
considerably larger than we might otherwise expect. 

Regarding immigration and emigration, it will suffice to note that 

in 1911 there were 350,429 immigrants into the United Kingdom 
and 454,527 emigrants from it. 

Further remarks on the constants will be found in Section X. 

(ii) W e ma} now proceed to consider numerical examples of 
proportional happenings. For this purpose it is convenient to adopt 
the suggestion made towards the end of Section VII (vi), and to take 
c as a function of an independent parameter y which is such that when 
y is positive KL is positive and a; always increases with the time, and 
when y is negative KL is negative and x always diminishes with the 
time (Section IX). We shall also write D for v — V and R for 
N -f- r. The principal constants of x and/may then be written 


c = D + ( Y + 1 ) R 


L 


Y 


Y + 1’ 


K = ( 


l 


Y 


1)R, KL — y R, 


Y 


(Y + 1 ) 2 


— 


(Y + 1) 

4y 
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If P„ is not 10.000, but some multiple or fraction of it, say 10,000 p, 
all we have to do to find t and i' is to add to or subtract from the 
figures given in the corresponding columns of the table the quantity 
1 oo f p yR. Thus, if P 0 — 1000, we subtract 2 * 30/y R. As y increases 
above 100. t and t' approximate more and more closely to the value 
loo, P 0 y] 1, t being always less than this and t' greater than it. 

The actual curves depend largely on the values of R. If we suppose 
tli e case of a human disease, of which the affected birth-rate N 
equals the natural birth-rate n, which is taken at 0*000083 daily, 
corresponding, according to the previous sub-section, to an annual 
birth-rate of 30 per mille , then R equals the former figure if r is so 
small as to be negligible. But if r is a daily rate corresponding to an 
annual reversion of 90 per cent., it will amount to as much as 0*0063, 
so that R — 0*0064. We will suppose that these two values of R 
are a minimum and a maximum. 

Now if we suppose also that c = 1, that is, that each affected 
individual daily infects or reinfects one other individual, then we 


hav 


e 




(1 — D)/R — 1. That is, if D = 0, 


Y 


12.000 if R is 


a minimum, and y = 155-5 (or roughly 150) if R is a maximum. 
Thus, if the infection rate is near unity, the value of y vail reach 
from three to five figures, and the epidemic will resemble the last 
examples given in the Table. 

In Example 1, if the population remains the same during the whole 
of this long period (which is, of course, never likely to occur), only 
about ten living persons ont of the 10,000 will be found unaffected 
at the end of it; and the ultimate and maximum daily number of 
new cases 1P„ reaches only 0-00083, that is, only one new case in 
3-3 years. 

If, however, R has its maximum value, there is one new case in 
about sixteen days. 

In Example 4, half the population is ultimately affected, but still 
/ has no maximum before its ultimate value. When, however, in 
Example 5, y increases above unity, this maximum begins to appear, 
and the ratio X between the maximum and ultimate values of / 
constantly increases with y. At the same time, t and t' diminish 
and approach each other. 
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If c is nearly unity, t' = 2-30 logy,, P 0 = 9-2 days, and/reaches 
its maximum in a short time, whatever the value of R. More 
generally, if I) = 0 . c -- (y + 1)R, and if y is large (say over 50), 

roushlv. and ? = He . log. P 0 . In other words, if y 


— 1 


is considerable, not only does L = 1 nearly, but the maximum 
(| c ) of / and the mode, or time(r') when that maximum is attained, 
depend roughly on c only, the population being fixed. 

If L > .! the f -curve slopes downward from its summit more slowlv 
than it rose towards that summit. Put x — \-\-\, then/ = c (-| — 5 2 ), 
and has the same value for equal and opposite values of 
But 





§- = |/-K(l-L)(l + 5 ), 


and is less when c is positive than when £ has the equal and opposite 
value. Hence, at two points of the /-curve whose ordinates 
are equal, the descending side is less steep than the ascending 
side. 


Since / 


c dx 
K dt 


c( 1 — L) x, and the first term on the right gives 


a perfectly symmetrical bell-shaped curve, we may call c(l — L) x 
the excess of the /-curve. When D = 0 and y is considerable, the 
excess varies roughly from |R at the centre of symmetry to R, when 
t is very large—that is, roughly, from |tl to l. Both of these become 
small when r, the reversion-rate, is small. 

The important result, therefore, follows that, if also c is near 
unity, and therefore y is large—that is, if the infection-rate is high 
while the reversion-rate is low—then the /-curve becomes a nearly 
symmetrical bell-shaped curve. -This, according to Dr. Brownlee 
(Section I), is just the kind of curve presented by epidemics of certain 
zymotic diseases. The prima facie inference is therefore that such 
epidemics may be n:holly or chiefly mere cases of proportional 
happening, as defined in Section VIl (i). It is even the case that, 
in such diseases, the fall of the curve is often more slow than its rise, 
as happens with the/-curve. 

When y is large and L = 1 nearly, the two changes of curvature 
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of the curve occur, according to equation (io. when 0.r 
that is, wheji (equation 53) 


f) 


\ 


*> 

O 


V It/ - loe, V„ log, (1 


,/ 


!)• 


- log P 


0 


1 • 32 and Ion. P 


ij i 


1 * U3 


tin; hit top being n can.*]’ to tin* sui limit than tin* former. 

(iii) Several methods may lie employed for calculating the curves 
of j\ /, / and F in detail. The simplest is to divide the period 
during which / varies Irom x {) to IJj into a niiinbpr (say 10) of equal 
parts and than to calculate the ordinate of the curve at each section. 
In equation 5S we used t to express this interval when x t} — 1/P 0 , 
but it may la* employed more generally to express the abscissa of 
tlu‘ centre of symmetry, when x n has anv value, so that 

It' t' a v * 

IvLt — log, (L;V 0 — 1). We now put t — rT/10. give to T the 
successive values h 2. 3. .... 10. and calculate the cor re. 5 
values of munelw ./ 1? ./‘ 2 , r 3 , .... and so on. Thus we have from 
equations 50 

J 'i — 



O’ 

o 


L{'x (L/> 


lu_ J’i — ^ 


l)"'- 1 } 


The root can be easily evaluated l>v means of logarithms, and the 

L <T 

value ot j; then obtained from a tabic of reciprocals (e.y. . as in Barlow’s 

1 allies). Owing to the symmetry of the -/'-curve, its values when T 

Vl— ^ * 

lies between 10 and 20 can be easily found—for example x Ui — L —x 7 . 
\\ hen u'e have ascertained the ordinates of x 7 we can quickly 
calculate those oty from equation 61 . and, unless the case-mortality 
is high, we can generally assume for short and sharp epidemics that 
the original population remains constant and that Z = xF 0 and 

^ It will be observed that, by this method, when t has been 

lirst calculated we need subsequently deal only with the values of 
L and of L f x ti — 1. Thus suppose that L 

1 — 10,000, and that x 0 = 1 / P 0 , then we have approximately, 

t - 0 
z « i 

- l 
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This is. of course, ail ultimate case where/is quite symmetrical— 
so that, for example/ 7 = cx-x i:( . 

Another useful method for calculating x and / is to transfer the 
origin to the centre of symmetry (equation oG). But if the ordinates 
are required for successive natural units of time, such as days or 
"weeks, m oidei to compare them with, statistics of epidemics m 

which such units have been used, then we must of course <dvc to / 

<”' 1 

the successive values 1, 2, 3, .... and calculate the ordinates direct I v 

f 

from equations 50. Here the symmetry of x will not he so usef 

■ 1/ 

because the centre of svmmctry is not likelv to coincide with anv one 

* v 

of the periods of time for wliich the ordinates have been found. It will, 
therefore, be generally easier to plot the curve by the method first 
indicated, to measure off the natural units of time aloimthe time-axis. 

4 






3 corresponding ordinates geometrically. 


Section IX 

(i) Hypomefnc Happening : KL Negative. —We have hitherto 
always supposed that the happening-element is constant during the 
whole of the considered period, but it is now necessary to examine 
the changes which will be produced in the functions by certain 
changes in this element. We here study onlv tin* cases of dependent 
happenings already discussed, and these only for hy pome trie 
happening. It is no longer very useful to take c as a function 
of y. 

In Section VII (vi). while considering the constants K. L. r. we 
saw from equations oQ that if KL is positive. ./• always increases and i> 
therefore always greater than its value when t —0. This implie> 
that c > 1) -(- R. When this mequabtv becomes an equality. I\L 
vanishes and d.r 0 jdi — {), that is, x remains near its original value 
when t is not very large, so that the happening may be said to he 
isometric [see (v) belowj. When, however, c < I.) p L* XL become- 
negative ; and it is seen from equations 50 that, as z~ wl ‘ then 
always increases with the time, x alwavs diminishes, and mint 


therefore be always less than x 0 . 


We call such cases hypometric happening. That is. we now suppo 


V 


e 
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that, after havin'* been greater than I) -j- U. c becomes less than it. or 

even zero. \\ hat, then, will be the nature of the curves ? 

I he solutions of equations 44. 45, 47. 51, 52 remain valid for all 

real values of e, and we have therefore onlv to interpret results 
aheadv known. 



It is first necessary to examine K, KL, and L more particularly as 

functions of c, but we shall do so only in the case when D is positive 

(injurious happenings) or zero. Then D. R are positive constants. 
We have 


K 


D, 


KL 


D - li, 


axis of c at an an<de 

of4o on the positive side of the origin. ° 


L = 1 — R/(c — D), 


or 


(L — 1) (c - D) 


R, 


L= w/V r * reCtangular h yP erI >ola, with asymptotes 

ardinatea. ’ d ^ “ the parallel to fche axes of co- 

The last portion of this L curve for r n l t? 

ircs; s ^ +r » Lr: r d :: 

^ ® * “oc® It becomes infinite when c — D, 


E 
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\vr haw to consider tin* two portions: 0 ^ c < I), tor winch [j is 
c and greater than ] : and 1) <c < I) R, for which 
\j is negative. 

(ii) First let c — 0 ; that is. suppose that the happening suddenly 
c(*as(*s altogether, while all the other elements remain the same. It 

i ' 

is advisable to use marked lett(*rs for the quantities after this event. 
W V take it th<‘n that when the original happening c ceases, f - - 0. 
and that the values wliicli ./\ Z. 1\ had reached up to this time are 
now (expressed bv ./*,/. Z 0 \ P () k while x . Zk denote these functions 
as (' increases. Similarly. K/ and L/ are now the values of K and L 
without tlio (dement c. Then, since K' = —I), equations 51 and 52 


beeonn 


P' = P 0 V r jl 


7 r _ 7 wv-ior 

/ j — /j 


1 \) 


L 


Ml 


„ - (IM P)/' 


) 


(si) 


c 



As V - P 


M f I - E - 


\ 


r. 7! decreases from its original 
alu<> Z,/ or Z. and ultimately approaches zero—unless the immigra¬ 
tion is greater than the sum of the other elements (as may happen in 
a war area)—and its diminution is more rapid when M. E. r are large. 
With regard to equation 81. we observe that P 0 V' is the natural 
increase of the population, living when the happening ceased, aftei 
that event: and the equation shows that P' equals this quantity 
less a function of I,’ which is zero when t' — 0, which increases with 
l', and which finally slowly approaches the limit .V/L . P„ - 
whim r is large. That is to say. the population which remained 
when the happening ceased continues to suffer from its effects lot 
some time afterwards, and then finally reaches the value 


I 


)' 


P„' s' r (l - .To'/L'). 


Here r,//J/ is always less than unity, for V is the value of r at the 
moment when the happening ceased, and this could not bait bun 
greater than L, which is less than unity [Section MI (vi) . win e. 
as shown in the previous sub-section, L' is greater than unity when 
c =, o. Hence P' is always less than and the difference 
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gives tlit? loss of population duo to the cases winch continue after 
the happening ceased. 

The fraction L7.r,/ occurs in the denominator of ,/ in e(juation 
o(). As it is greater than unity, the. coefficient of 2 K 1 ' is positive, 
so that x' always diminishes as t' increases—-as could hav<* been 
ill ready inferred from equations 81 and 82. 

(iii) The next case is when c' is some quantity between zero and 1) 
(which is generally very small). For tin* study of this case, the 


first of equations 50 may he written 




' I x o' 

0 w 


(1 


x 


0 Y_ — K'l.'/' 


-1 


IV 


(8:i) 


which diminishes indefinitely as /' increases. If c ‘ = 1) ami L' 

is consequently infinite, this becomes .r oS K ' ,r . which vanishes. 

The values of P and Z' are easily obtained from equations 51 and 

52, if we remember that the indices D/K' and c'/K' are now negative. 

When c' = D, the expressions for P' and Z' become indeterminate, 
but can be easily evaluated, and we have 


P' 


T> 

x 0 - 


where q = x<) 'D (l 


- ft/' 


)/R. 


(iv) When c' lies between D and D -|- R, K' is positive and L' is 

Oivl £ __ r __ 


negative, and varies from 


to zero. Then it will be seen from 


-->vm ut* stMMI irom 

7 Wt, . # f 83 * h f *' “ stiiI ro si ‘ ; ™. »nd, « in the previous eases 
unmmshes imiofimtcly as f increases. For it may be written 


X 


xJ le~ K ' Lr — X ° /- — K'L'T' i 1 

1 L' ' ~ 1 ); ’ 



where K'L' is negative. For P' we have 


x, 


P' = P„'e ir <|l — (i 


K'j/r 


) 


-D/K' 


(85) 


which is also always positive since L' is neoativp n 0 ■ 

TTnT when r is p ' - less , 

aftra th° ti * ° rlgnlal l* a PPenmg continue to be felt by the population 
after the happemng itself has been reduced from c to o'. 
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(v) Y\ hen c' = D + K exactly, L' and KL' — 0, and, when t* 
is small. :r — .r,/ nearly, and the happening may be said to lie 
isnntcln<\ Tlie values of x and V' become indeterminate, but can 
be ascertained by finding the limits wlien 1/ = 0 of the expressions 
in equations 84 and 85, so that 


r f 


J 0 


1 + rAKY : 


P' = P 0 V r (1 -|- .r ( /KY)- 1,/, \ 


Or the same equations can be obtained by integrating equation 
47 after putting L = 0 in it. and then integrating both sides of 
equation 44. In this case also, then, x' tends to decrease as t' 
increases, and finally approaches zero, when t' is very large. 

(vi) If K'L' is positive, however small it may be, x‘ always in¬ 
creases, as already seen in Section VII (iii). 

It must be remembered that, though x , the proportion of affected 
individuals, diminishes, this does not mean that new cases cease to 
occur. On the contrary, unless c is absolute zero, /' is the same 
function of x as before. And rx' expresses, as before, the reversions 

of the old cases. 

"When D = 0, the phase described in (iii) above is suppressed, and 
we begin, when c = 0, at the point when K/ = 0, IvL - R? 
and L' = — « , and the functions will be easily understood. 


Section X 

(i) Parameter Analysis .—We have now considered x. f, and I as 
functions of the time, but it remains to examine how they vary if 
we give different values to the parameters h, ?, v, V, n, m, i. e, N, II. 
I, E. The symbol l has been used to denote the limit of / when t 
is very large. This is an important quantity in disease happenings 
because it gives the number of new cases which continue to occui 
when x has reached its limit, that is, after epidemic manifestations 
have ceased. It is, by hypothesis, these final new cases which 
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v> 

j * ► 


keep the infection alive permanently in the population ; we may call 
/ the endemic ratio. First consider l as a function of r 


cl 


< r 


- (21 

Iv 1 


1 ) 


0 / 
K- (c 


D — 2R) 


I inis / is positive when r — 0 (if IvL y> 0) and increases with r until 
this reaches a value which makes L = and l = \c [Section YII (vi) h 
-Viter this l diminishes as r increases. But if c — D > 2R then r 


never reaches this value and l always increases. Thus / considered as 

a function of r has forms similar to those of/considered as a function 
of t. 


1 he behaviour of l as a function of e is of the same nature. A. 
similar procedure applies to/and P. 

fhe effect upon the total population P due to variations in the 
parameters is important. The fundamental equation is 44 : 

dPjdt = vP — Dap. 

-Vs P and x are always positive, the sign of dPjdt depends on v and 1). 

If D = 0. the happening is equivariant and the natural change of 

population due to v is not affected. If I) is negative, the happening 

is beneficial and the natural change of pojmlation is augmented by it 

If D is positive, the happening is injurious and the natural change 

is let need. If v — Dai is negative, especially when x = L, the total 
population null diminish indefinitely. 

It is found that ?P/c> is always positive ; that is. in infectious 

is hZfi 3 , qlU ^ W l0SS ° f afleCte<,ness - including loss of immunity, 

- v ? a ' T “ s ™ a - v seem surprising until tve reflect that, as 
, e longer the affectedness continues in the individual the 
greatei mil be the loss of life due to it, since we generally suppose 

m such cases that M > „ and that these elements act during noj 

Of iol^to^Vr 2' thC affeet n ^ ° f ^ “ i0 ”- 

affectedness,^ including in^umly^ S 7 °“ ^ 
e n,//-:1 s Zei '° , negativo " and the higher the infection rate 


the more injurious is the happening. 
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f r ' ,s zero wJ,e11 * = 0. It generally increases at first with f 


and afterwards decreases 


it is negative when x 


} 


= L provided 

r 0 14, h- In other words, in these cases an 



increase of \ is beneficial at the beginning of the outbreak, 
injuiious later 012 . 4<ov r Y = h — M -f- I — E and diminishes if 
Ibc case-mortality 41 increases; in fact rP fdU = — rP/aV. Thus 


au increase m the rase-mortality is injurious at first, but may be 

i' ht ^ ml, because an affected individual who dies ceases 


to be infective. 


oo 


PAliT Til 


Section XI 

(i) Variable Happening. Tn Section IX. we commenced by 
touching upon the necessity of studying the effect of changes in 
the happening-element, but there dealt only with the case of hvpo- 

nietric happening. AVe now proceed to examine other changes in 
this element. 

It is unlikely that in any infectious disease; the iniectivitv-element 

t 

c will always remain constant. It may remain practical] v such for 
some length of time.; and the hypothetical results of this condition 
have therefore been analysed above in detail : but nature abhors 
a straight line, and we may infer from general experience that 
changes are sure to occur from time to time. 

Section I. changes in infectivity mav be due 


As suggested in 


(a) to the action, so to speak, of the infecting organisms themselves 
--winch is difficult to believe in because it involves the conception 
that organisms living in large numbers of different hosts should, as it 
were, make a land of concerted and simultaneous effort ; and 

(b) to changes in season, climate, or other environment which mav 
(i) increase or diminish the infective strength of the infectiim 

Sranothe r r . <2) “““*** " 

l " lt ,ileSent . we .°“ u P roc<!ocl only on the supposition that the change 
has occurred, rvrthout attempting to enquire how or why it Iv, 

■ “ es :r r “■ 

lartre Tf + i,. i . . iauals 1S z o 5 which may be small or 

the latter figure- bnt if m i J^ais, Z, wffl increase towards 

° ? e happening is isometr ic (KL = 0), 
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is iiy|)f,metric- (KL < «). Z will diminish to the* limits discussed in 
Section IX. \\ hat will happen d the infectivitv is suddenly cliimuei 

4 ( 1 

during one unit of time from c to r', and thereafter remains at 
th(‘ latter figure { 

4 ~ 

Id v > c\ the ./-curve suddenly becomes steeper : itself is con¬ 

tinuous. but its gradient is discontinuous. The /-curvi 1 is djscou- 
tinuous. and has a sudden rise at the moment when c changes ; 
thereafter both curves j)i‘oceed according to tlie e(jnations with d 

substituted tor (\ It at the moment of discontinuity more than half 

* 

the population is affected, so that /is descending, there is an epidemic 
manifestation, a sudden rise followed by a decline. Both L and / 

4 . 

are increased. 

These are probably just the conditions which applv to manv 
zymotic diseases, in which most of the population appears to remain 
affected, that is immune, for years or permanently. In those cases, 
at tht* moment of discontinuity x is nearly L. and L is nearly unitv. 

4- %/ ' L t 

and the augmented infectivity d acts merely upon 11 re small residue 

of non-afiected individuals, let us say, newly-born infants,, or children, 

or others who have hitherto escaped infection by accident or by 

residence in remote villages. It is well known that at the beginning 

of the war recruits from many highland villages we re attacked 

bv very virulent measles on joining camps in the lowlands. On the 

other hand, the same sort of epidemic may occur in numbers of 

villages and small towns in which the considered infection has died 

out from want of material, and so allowed of a gradual accumulation 

of fresh material, which is then suddenly fired (so to speak) by the 

entry of one or two infective cases from outside. 

* 

It is extremely improbable, of course, that the infectivity actually 
makes a single per saltum change : but we may still apply the ccpia- 
tions already obtained if we give to c a series of progressive dis¬ 
continuous augmentations for successive units of time, followed by 
similar discontinuous diminutions until e returns to its original 
value. In this case the /-curve will again show epidemic manifesta- 
lions (if x > f). and will then decline again towards a. lower value. 
This is one way of approximating to a continuous change in <\ 

(iii) Continuous Variation.—We may, however, obtain new 
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equations on the hypothesis that the inl’ectivitv r is a continuous 

1 

function ol the time. Let us suppose that r — </>' (/). wIhmv 

<y (0 - <t<f> (/) '/L Then in thr fundamental equation 47. namelw 
(hr di -- K.r (L — ./ ). we have 

X — <f>‘ (/) — I) ami KL — <f>' (/) - f) — 


Put . 


1 //. so that dt/ dx 


K \j)j .- K : then 


// ( ‘-\p 


K \jilf 


( K exp 1 IvL dt i dt -- constant 


by the well-known solution. As however. K 


KL l( 


K exp |i KL dt I dt 


■i ' 

exp KL dt —i- I* 1 1 (>.\p 1\ \j(h MU ; 


h 


ence 


a 


] -1- exp ( 


KL dt i It 


• o 


f 


O 


Tt 


< exp i K L dt l dt - - if 

* 0 


1 


4) 


Tl 11 is negligible. as may often happen in short and sharp epidemics 
with long immunity, we may write, for many diseases at least. 


KL dt = <f>{t) — D/\ 


x 


1 


1 + const 




„ - Ain -I it • 


C 


We have at present no knowledge of the nature of the function 
= <f> (t) or eyen if there is such a function. But if then* is. we 
may suppose that, for positive values of t, <f>' (t) is small when / is 
small, increases at first as t increases, reaches a maximum for a 
certain value of t. and then declines to a small value again when 
Os larger—as is observed in many phenomena. In other words we 
should expect a symmetrical or non-svmmetrical bell-shaped curve 
such as those with which we are familiar in statistical work 

(iv) Graphical Treatment of x .—We may also suppose that c 
is given in terms of the time t, not by an analytic function, but bv 

a curve whose equation is unknown. We can then determine the 
natm e of x and f by graphical methods. 

Since djpit = Kr(L — a), there could only be a turning value of 
^ when K = 0 or x = 0 or x = L. But when K = 0, then KL ^ — B 

T ‘'° es vanisl ' : » = « the happening ceases 

Itogethei , the only turning value of x is when x = L. 
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^ tnv ** llce O' definition x lies between 0 and 1, we can only have 
— L if L also lies between 0 and 1 ; the diagram of IX (i) shows 
that this only occurs when c > D -f R. Hence unless the given 
r-curve rises above the level D -j- R. the corresponding a-curve has 
no fluctuations, and either rises steadily to 1 as its limit, or sinks 
to 0. But wherever the c-curvc has a portion rising above I) -j- R. 
we can ]>lot the values of B = 1 — R/(c — D) corresponding to this 


portion, which form a cap. part of the L-curve. lying entirely between 

’ r v ' 1/ 

the levels o and .1. If the e-curve fluctuates, there are a series of these 

L-caps. which are the only parts of the L-curve that the .r-curve 
can possibly cross, 

^ow. whenever the ./-curve meets an L-eap, x — L. dr (If o, 
and x has a turning value. If the ./-curve descends to meet the L-cap 
from outside, it begins to ascend as soon as it is inside ; and goes 
on ascending till it meets the opposite side of the same cap, when it 
passes outside and immediately begins to descend, until it meets 



Fig. 1. 
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oU 


another L * up Till- /-curw descends steadilv except where if is 
iiUMile -« raj# and a- -ron the L-caf#^ are plotted, the general 
1 hmetef of the curve is fairly apparent* If c has a definite limit 
' , when I x and this limit is less than D -f» R, then the last 

1* * f' 1 * « (r»«» d , is ultimately descending and tends to o. But if 

«r O 


/ 


' D R, then L has a limit h t and x tends to the same limit 


il ‘ L and may approach it either from above <»i from belowx 


D), t Is >r K (r — Dj-, and is positive 


Since L = I — H/(c — 

I, i j ~e.> amj fall > with f and the top of an L-cap corresponds exact I y to 
a maximum of r. But a maximum of x occurs where it crosses a 
descending branch of L : hence the maxima of x occur somew hat 
hlter than the corresponding maxima of c. 

Suppo-#* foi example, that r is periodic. rising and falling regularlv 

ever} jrw, and cro->mg the level It j- R twice a v- r. TIiqd L 

• » 

lias a ncri> ^ ..r i-ups with *<pial <<pii<Iistant maxima of ln-i^ht >., saw 

>ruf Irom a minimum of < Iw-tw n two L-caps, ami suppo'i- tin; 

ilu* - <„ of x to be given. 

llit-n / bfi'ins to «I*->o.ik|. If x u is lai "*• i may pass over the top 

{ v a ■ 

nf one or more LrOaps, dfiaoending all the while, but, sooner or later, 
ir niet-r > an L-cap, / - * - ^ m-id*- it and begins to rise. It thereafter 

meets every etp, for it leaves each cap at a height less than X, and, 
race it deman ds between the caps, it cannot pass over the top of the 
uex Kv«rv v ear x na^ and falls, the maximum occurring later 

than that of i or JL, and we have a yearly epidemic. But there is 

no r* a ‘»n to mulo- these epidemics equal; though, if any tw r o are 
qua! tli* > all mist be. For if / has equal values i*»i corresponding 
1 s of the periods, the whole rourse following these points is the 
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same ; in other words, if two periods are superimposed, and if the 
two parts of the x-curve have one point in common, they coincide, 
but, if x meets any L-cap at a point lower than the preceding cap, its 
whole course inside that cap is lower, and lower still in the next 
cap, and so on, tending to a limiting form which is periodic, and which 
may he the axis itself, the happening dying out. But, if x meets any 
L-cap higher than the last, it meets the next higher still, and so on, 
and tends to a periodic limiting form of finite fluctuations. 

In order to plot x more closely, we have to replace the fundamental 
differential equation by the approximate difference equation 



where K, L are calculated from the value of c measured off the given 
curve, and x 0 is given ; or by some modification of the method. 

(v) Graphical Treatment off.—U c is given and x has been obtained, 
we have next to deal with /, which can be directly calculated from 

the equation 


/ = cx (1 — x) = c {£ — (x — 4) 2 } = c {£ — (J — *) 2 }- 

Hence/ always increases with c : as x increases, / increases if x 
and decreases if x > If we wish to know the general behaviour of 
/, the following Table is useful; it allows for a single bell of each of the 
c- and x-curves, the maximum of x occurring later than that of r. 



i 

i 

h 


/. 

c. 

* 

if * < J. 

f 

if J* = b 

t 

minimum 

decreasing 

decreasing 

minimum 

! 

increasing 

decreasing 

• 

1 

increasing 

L 

increasing 

i 

minimum 

j increasing 

i 

! increasing ; 

; l 

9 

increasing 

increasing 

i 

increasing 

increasing 1 

# 

maximum 

increasing 

1 

i 

increasing 

P 

, maximum 

decreasing 

increasing 

L 

• 

decreasing 

decreasing 

| maximum 

decreasing 

} 

decreasing 

decreasing 

decreasing 

i 

decreasing 

decreasing 

* 

minimum 

j decreasing 

decreasing 

minimum 


if -** > 4* 


increasing 

increasing 

increasing 

% 

decreasing 

decreasing 

decreasing 

* 

* 

increasing 
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In using this, we follow/down a column until x passes through t he 
value J, when we cross to the next /-column. 

Since in each column f is increasing in the third row and decreasin'' 

i / t ' ^ 

in the seventh, there is bound to be a maximum of / between 
them, which occurs before the maximum of r if x > J, and after 

if./: < t. 

Recrudescence .—It is seen that if c: and x are both decreasin'' and 

- we cannot say what j is doing, as the two causes of its change 

act in opposite wavs. It is possible for / to have a minimum 

followed by a short rise and a maximum, and. if ./■ then becomes 

equal to A, / falls again to another minimum. Tins would happen. 

1 o* example, if c suddenly became constant after decreasing, before 

x drops to h Then we have a slight recrudescence of the epidemic 
near its close. 



Since 


<n 

dt 


( 


f dc 


Pc dt ?x dt 


f dj 


and in the case considered dc/dt and dxjdt are both negative, 


dt 


x (1 - x) 


dc 

dt 


+ c (2x - 1) Ivr (L ~ ./•) 


and is positive, so that/ rises again, provided dc/dt is small enough, 
but, as ,r decreases, the factor 2x — 1 decreases'and changes sigm ; 
t ie positive second term in df/dt becomes negative, and/ sinks. 

In fig. 3, c rises and falls in a symmetrical bell and then remains 

constant; L does the same ; a- has an unsymmetrical bell, risim. 

more steeply than it falls, and sinking to the same limit as L • and 

/shows an epidemic with slight double recrudescence before sinking 
to its limit. ° 

Tke graphs of c. x, and L are on the same scale, and that of f on 
lu times the scale. J 

(VI) / Assigned .—The last sub-section started by assumino- that 
C IS given in some experimental way. But c is not given except by 
conjecture and the method is useful to test whether any partLlar 
conjecture ,eads to numeriea! results which agree with "he“ s 




















n <IS s 



r If.'. -J 
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of a< ?tial t*pidemirs. Now tin* most extensive and ac<-ui*atc records 
an* the nofitications of fresh casr.s and of deaths mad** bv the medical 
oliic'ojs ot health : in particular, the records of measles in Aberdeen 
for tie* ifo vears Is Hll l!M»2 are wrv suitable material. thour'll even 

* If ^ 

t here the numbers are not larjie enough to be quite satisfactory. 

I he i [eatIi-rates do not correspond to anv of the symbols we have 
i nt rod need so far ; tor M / is the number of deat hs ot affected persons 
per unit time. a different tiling from ttie number of deatlis due to the 
affection. W e return to this in Section Xf I (i). 

\\ e start therefor.- with tin* fresh assumption that the notifications, 
that is. th*- values of F. are known. Suppose also that I* is given, 
either by interpolation from census returns, or because in a short and 
mild epidemic I* is nearly constant. Then /is known as a function 
o! /. The fundamental equations are best taken in the forms 


h dt ~- f - (|i -y K — Da), 


c = f! x (i — x) 


from whieli a and r can respectively be determined, if reasonable 

•*>tunates can I>e made of the constants JJ, R. 

When a, / are found and plotted, we can construct r graphically. 
I* is the point (f.J). Along the vertical axis mark off 


< >A 


x. 


OH — l o(< 


1 


r- * 



AH - l 


x. 


AC 


1 


le t 1’A meet the axis XO in V, : let \'.H meet PX i 


\ and let V.,( ’ meet f J X in |{. 


in () anrl XA in 


Tlua 

and 




HX/AC = QX/AB, 


QN = fix, 



f J ' (1 — x) = c . 


S ° I]'** H on the c - curve corresponding to P on the /-curve. 

If D and U can both be neglected, then 


rt 


* o 


F dt 


z 0 . 


In this case ar steadily increases whatever is the law of / and the 
area under the /-curve alwav* ^ i J ’ the 
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In a iiiv.it manv cases, the /-curve is a slightly unsymmetrical 

* * 

hrll. l isini* move steeply than it falls. Isow 


Ml 


•'• 2 ) v = %- x l 1 “ a ‘) —f *( l 



dt dt 


ami df dl is positive for the first, smaller part of the range then 
vanishes and is negative for the rest of the time. The second term 

is negative if jc < and positive iix> 

Now if /„ is small, at the summit where df[dt — 0, x < 2 an 

Tn this case the c-curve reaches its summit 

*■ t^tlrnTthat of /, and is more unsymmetrical. But if «. is 
Zlr 1 - Sr Changes its sign earlier, and the summit of the c-curve 
i„,er in consequence. These results are an extension of th 

same case where D, R, are neglected, /is perfectly symm 
" Hu,"do not wish entirely to exclude reversions, we assume 


D = 0, 


R > 0. 


Then 


dx/dt -f- R® —/’ 


x 


-1« 


j l ft Kt dt+X 0 e- lu . 


. a ,(T — t) which is the simplest function that 

rise s and falls, vanishing at f = 0 and t = T, with 

l = l T. Then 


. - fi{R (T — 0 + 2} ^ S (RT + 2) (1 

Or we may allow lor asymmetry by an 
f= at (T — t) s'M which still vanishes at t 

has a maximum before t = \ T ; 


- E ' E ‘) + *o e_Ri - 

exponential factor, 

= 0 and t = T, but 


at 


- Tn 


(R - M 


{(R — X) (T — t) 4- 2} 


—At 


S-J, «B + X, T + 2) (.- 


-1U 


) + Xq 


—Bt 


X 
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If / is periodic, let 
/ = a + b sin pt, 


x 


a 


R 


7(1 




) 


where a — b 


b 


P 


TT 2 ( r sin pt — p COS pt + pz-' u ) -r x 0 z 


_ - lit 


In all these, the negative exponentials decrease rapidly ; in the 
last, x oscillates about ci/lR. as mean, the phase lagging behind that 

°1 j by — l an _1 ^ • Then c has the same period. 

(vu) More General Equations. -In many diseases, for example, 
those m which there is a period of incubation, the infectivity is 
variable, and depends not only on the absolute time but also on 
the duration of the case , which is the time elapsed since the case was 
first infected. In order to take this into account, c must be con¬ 
sidered as a function of two independent variables, t the time, and 

S * h ® duration ; fluctuations of the infectivity of the disease’as a 
w ole cause c to vary as a function of t ; episodes in the course of a 
single case cause c to vary as a function of », very possible a dis¬ 
continuous function. For example, an incubation period of , 
days is represented by taking c = 0 when a < q. 

delencfonTanf 116 ?'I 7 ’ SUPP ° Se ^ V ’ h E ’ M > N ' c - ^ D > K . 
depend on t and s ; but v, i, e, m, n on t only. 

t and f+" l be ^ " Umber ° f £reSl1 04863 ooc “ r ™g between times 

Let F t+S s dt be the munber of the batch nf F /it , . , 

survive, having neither died emitted ° ' “ \ a "' hl ° h 

t + s, bemg then of duration s. Thus F i« +h. i \ * e 

existing at time t whose duration is 5 Th t nUmber of cases 
existing at time t is • e total number of cases 


Z 


n f, 

J 8 =* 0 


ds. 


( 86 ) 


As before, P — Z + A, and the number of fresh 


cases is now 


F 


t, 0 


A 

P . a = o 


cF,, t ds. 
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The rate of change of the particular batch of cases which started 
at time t 0 is given by 


d 

dt 


F,.,.,. = (V — N — r) F,,-„ 


which may be written 


where 


iL _j_iLw 


= (V — R) F, s , 

— t t 0 . 

The fundamental equations now take the forms 


dP 

dt 

dZ 

dt 

dA 

dt 


r oo 

vA + VF, s ds 

* 8 = IJ 


OO 


J = 0 


(V R) lh. .<? ds h 


0 


vA 


CO 


. -= 0 


HF t . t ds-F,, 0 J 


(87) 


( 88 ) 


These equations are not independent. If 87 is integrated with 
respect to s from 0 to oo , since F tt w — 0 we have 



rao 






00 


J S — 0 


(V - R) F'.,ds, 


which is equivalent to 88 by virtue of 86. 

These equations as they stand are quite unwieldy, but they are 
capable of being applied to a large number of simpler cases. 

If V, R, are independent of t and constant for a particular range 
of $, say s — qi to q 2 , then for any intermediate value of s. from 87, 

F — F _<v-it) (*-?,> 

t, $ * t—q.- q l w ? 

and if V, R have other constant values Y'. R' from -s = q 2 to q 2 . then 
for any value of s in the latter range. 

F = F _(v'-K')«-«,) 

- 1 L 8 J t—q*. q* 

_ ir e .(V'-K')r+<V-V # -U + lt') 7 ,-.<Y-J!) 7l 

- 1 t— ?i- ?i w 

and so on. 

The simplest case is when c is constant from s = to q 2 - zero 
for all other values of s, and independent of t, while all the other 
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tion is a constant M; part of this, m say, is due to other causes, and 
M — m is due to the disease. If m is constant, so is M — m, and 
the deaths from disease are simply proportional to Z and their ratio 
to the total population P is proportional to x. It therefore increases 
steadily ; this does not agree with the records. 

It is more reasonable to assume that the deaths from the disease 
all occur soon after infection and that their number is proportional 
to F ; let it be |aF. Then its ratio to the total population is y.f, 
and rises and falls in proportion to /. This is much more nearly 
true. But now the fundamental equations are altered. We assume 
that the (aF deaths are in addition to MZ deaths due to other causes, 
and that they all happen before the cases have had time to infect 
an appreciable number of other people. Then 

dP/dt = VZ + vA - (aF, 

dZfdt = (Y - N — r) Z + (1 — |a)F, 

F = cZA/P. 

As before write 

Z = xP, A = (I — x) P, F = cx (1 — x) P, 

v — V — D, N -j- t = R. 

The result of elimination is 

dxldt = - {D + R - (1 - n) c} £ 

+ (D — (1 — 2*j i) c } x 8 — (acx* 

— — jj .cx (x — a) (x — (3), say, where a. > (3. 
of which the integral is 

g -MCaP (a-0)« __ /£<*-0 (x - a) 3 (X - (3) -a /{X 0 a ^ {X<> *) ^ 

The form of t his curve depends on the signs of a, 13- 

Let Q = |AC (x — a) (x — (3) 

= [xcx 2 — {D — (1 — 2 (a) c} * + D "h ~^ ^ 

, • - /, \ / P 4- R, then 

Case (i).—If the last term is positive, (1 — (a) c < ^ t 

«, (3, have the same sign ; and if x, is a small positive quantity, 
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than 


or 


- , - , — | J, dx 0 /dt is negative, x begins to decrease, Q cannot 

change sign and x decreases steadily to zero as its limit. 

I* the last term in Q is negative, (1 — jr) c 

When x = 1, Q = R and is 


Case (ii). 


> D -j- R, 


then a is positive and |3 is negative. 

positive ; therefore 0 < * < 1. Now dx 0 /dt is positive, x begins to 
increase, but it does not reach « till t = oc : " 

j 

sign and x increases steadily to a as its limit. 

"Ik 'T" /> 


hence Q cannot change 


0 

© 


Now/ = cx( 1 — x), and as x varies from 0 to 1, / varies from 0 to 0 
in a bell-shaped curve with a maximum at x = -t. 

In case (i) / decreases steadily from/ 0 to 0, in “a curve of the same 

general character as the avcurve, for the factor (1 — x) does not 
vary much. 

In case (n) if a < -J, / increases steadily to cx (1 — a ) as its limit 
utit a > -t,/has a maximum \c when a; = \ and then sinks to its 

limit COL (1 — a). 

= , D ’,? aT A h ° th Ue * IigibIe - Q = c (a - 1) (^ + 1 _ ^ and 

a .. ’ “ en / ^ as a maximum and sinks to 0 as its limit. 

(ii) Second Modification.—A popular assumption about infectious 
c iseases is that the case ceases to be infective after a certain time 

;“° tWe ° eaSeS *» be *“«« d-th from the disease, bui 
cma ns immune permanently. To represent this, suppose the 

population P to consist of three classes : P=A + X + Y xh 

A are unafiected. X infective and immune / LLIi 

simphcity, nedect all “ i ^ . mmune. Tor 

MX deaths from the disease 8 per^ afe 

' = CXA/P ^ ^ame^r^ 


dP jdt 


MX, 


dY/dt = rX. 



dA !dt 


cXA/P — (M + r) X, 
cXA/P. 


Hence P decreases and V increases stead,ly. SuppoM Y< 


0 . 


dY/dP 


r/M, 


dA/dP = C /M . A/P, 
(l/M)dP dt ~ X = (i _j_ p 


Y = (^/M) (P 0 — P) ; 


A 


A 0 (P/P o) mi 


t 


MJ 


(WM) P 0 —A 0 (P/P 0 )"'I ; 


. . rfl> /« > + r/M) p _ ( r/JI) p o _ Ao (p/ p t ylx 
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which is integrate in finite form if C /M is rational, and leads to 
algebraic functions and logarithms. 

In the final steady state, X = 0 and P is given by 

(1 + r/M) P—(//M)P 0 — A 0 (P/P 0 ) elil = 0, 

which has a real root P„ between P 0 and 0. A definite number of 
persons A„ = A 0 (P x /P 0 ) s/M escape the disease altogether. 

Now, dX/dt is of the same sign as cA — (M -f- r)P, and is positive 
at first, provided c > (M + r) P 0 /A 0 ; in this case X, the total 
number of infective persons, rises to a maximum before sinking to 0. 

dx/dt is found to be of the same sign as (c — M)A — rP 0 , which is 
positive at first if c > M + ^P 0 /A 0 , and steadily decreases. 

dF/dt is of the same sign as cA — (M -f r)P — (c — M)X, and, 
if c > M, this is less than dX/dt; the maximum of F (if it has one) 
is therefore earlier than that of X or of MX, that is to say, the 
greatest number of deaths occurs later than the height of the 
epidemic, which is the time of greatest number of notifications. 

For example, let 


r = 0*01, M = 0 *02, c — 0 • 04 , x 0 = 0 • 0102. 
Then A 0 /P 0 = 0*9898. Let P/P 0 = p, 


M/ 


CP 


dp 


,-i 0*9898p 1 2 — 1*5 p + 0*5 


0*52 


log. 


V 


1 + 51s 


— 0 0104 / 


2 • 02 4- 49 • 98e 


- 0010 « 


( 2 • 02 p — 1 — 0*02 1 

[0 *98p — 1 1*02 J 


P„ = £P, A w = ^P* approximately. 


X is stationary when A = fP or 1/p = 1*32, leading to t = 300 
approximately, the time which elapses before the maximum of the 

mortality. 


Section XIII 

(i) Discussion.—In the foregoing work, very little attention has 
been paid to the actual values of the constants. We have rather 
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been preparing apparatus which can be used in a great varietv of 
special cases. 

It will be seen that the number of dependent variables is large. 

and that a fluctuation in any one can be accounted for in a meat 

many different ways by supposing that there is a suitable adjusted 

variation in almost any one of the parameters. And the variation 

of any one parameter when it is over a range small compared with 

the mean value can be represented approximately bv one of several 

analytic functions, which would lead to the most diverse results 

if the variation were to continue over a wider range. It is therefore 

extremely rash to rule out a priori any one possible cause of fluctua¬ 
tion of any one phenomenon. 

The cases already considered have led exactly to the series of 
curves required by the facts : 

1. The steadily rising curve of a happening 
permeates the whole population [VII (hi)]. 

2. The symmetrical bell-chaped curve of an epidemic that dies 
away entirely [VII (v)]. 

3. The unsymmetrical bell of a new happening that begins with 
an epidemic, and settles down to a steady endemic level [VII (v)l 

4 The periodic curve with regular rise and fall due to seasonal 
disturbances fXI (iv)]. 

5. The more irregular curve where there is recrudescence before 

ho end of an epidemic, or where outbreaks differing in violence 
occur at unequal intervals [XI (v)]. 

rise anYfflfiT ’ “ WaS Stated in the Production to Part I, the 
rise and fall of epidemics as far as we see at present can be exnla ned 

y T , 6 ge "f a ! a "j s of happenings, as studied in this paper. 
rem“ ° ' W ^ " Sed ° al1 ' for the follow* 




o 


fT“°" VI1 ' «“ case-mortality and 
that this is strictly true “nd it t “ iS ^ U '“‘ teI 3' 

values are highest at som A « * “ ? UCl1 m ° re Probable that their 
then gradually die away. ‘ Peri ° d afte1 ' infeotion and 

B. In Section XI some attempt is made to meet this objection by 
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introducing the case-duration as a second independent variable, in 
addition to the time, but there are innumerable other considerations 
that may also be taken into account, such as the age, sex, and social 
position of the patients, climate, temperature, and, in the case of 

insect-borne diseases, the habits of the carrier. 

C. Except in XII (ii), the whole population has been divided into 

two classes only, the unaffected, and the affected who are considered 
to he both infective and also immune. But it is again unlikely 
that this is strictly true, and, in order to represent the facts 
accurately, a much larger number of classes would have to he 

considered. 

D. Although it would probably be easy to write down equations 
on any given hypothesis as to the number and behaviour of the 
variables, these equations, like those of XI (vii), would probably be 
quite unsuited to numerical work, and we should have to fall back 
on laborious approximate methods such as expansion in series or 

mechanical quadratures. 

E. The quantities x and L have been treated as continuous 
functions of the time, but, as they represent numbers of persons, 
they can only take integral values, and must be discontinuous 
unless they are constant ; in the same way, / and l can only van 
by integral multiples of 1/P. This is a serious drawback when the 
numbers are not large. In a small village, the average number o 
cases in the period between two epidemics may be only two or three, 
and, if these are removed by any chance, the whole course of events 
is entirely changed ; there is no one to keep alive the infection, and 
a new generation may grow up entirely non-immune. Then the 
entry of a fresh source of infection has the same efEect as if c rose 
suddenly from 0 to its normal value. In fact, in the case of sma 
populations, other than the most probable values may often obtain, 
so that the whole of this book is applicable to large numbers rather 

than to small ones. 
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